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Abstract 


This report is concerned with estimating the high- re solution amplitude and phase 
parameters of the spatial Fourier transform of an object illuminated with incoherent 
radiation and viewed through a turbulent atmosphere. 

The Cramer-Rao technique is employed/to lower-bound the variance of any unbiased 
estimator, and receiver structures that approach these bounds are then exhibited. It is 
found that significant improvement over many existing systems is often possible, pro- 
vided one is willing to pay for increased system complexity. The reduction of the 
processors to somewhat less efficient but more practically attractive forms is then 
considered. 

Some experimental results are presented which illustrate the effects of turbulence 
on wave -front coherence. 
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L INTRODUCTION 


1 . 1 STATEMENT OF TEE PROBLEM 

It has long been recognized in astronomy that angular resolution is controlled 
primarily by the refractive index inhomogeneities in the atmospheric path. These 
inhomogeneities result from the' continuous turbulent mixing of air from layers of 
differing temperature. 

A wave that propagates through the atmosphere suffers random fluctuations in phase 
and amplitude. Thus the image formed by a telescope is diffused, thereby reducing the 
resolution of detail on the object. This problem, the bane of astronomers, is often 
referred to as 11 poor seeing". A familiar example of this phenomenon is the blurring 
and dancing of an object viewed through an optical path containing a source of strong 
local heating such as a radiator or roadway. 

Even upon the best mountaintop locations, it is doubtful that wavefront coherence is 
maintained over a large enough area and for a long enough time for the full diffraction 
limited capabilities of large telescopes to be attained. These large apertures are of 
value mainly for energy gathering. 

In many situations such as satellite surveillance one is even more constrained in 
the choice of observing location and time. Average apparent spreading of a point source 
is often from I to 6 seconds of arc (1 second of arc = 4. 85 |Jtrad). Daytime observa- 
tions are usually the most severely affected. 

In recent years, there has been a strong revival of interest in atmospheric optics, 
which is partly due to the development of laser communication systems. The resultant 
cross -fertilization of physical optics, system, and communication theory has given rise 
to the application of signal analysis and processing techniques to the imaging problem. 
This was the subject of a National Academy of Sciences study in the summer of 1966 at 
Woods Hole Oceanographic Institute, which stimulated the present research.* 

Much of the contemporary work has dealt with improving the quality of images formed 
conventionally with a telescope and photographic recording. In previous applications, 
however, statistical communication theory has made its most significant impact when 
released from the narrow confines of a specified detection scheme, and allowed to seek 
instead the basic limitations of the problem within which all schemes must operate. Thus 
through such techniques as maximum-likelihood estimation, the best receiving structures 
can be found, and at times may be fundamentally different from the hitherto employed 
techniques. Also, once the basic performance constraints have been found, we are able 
to more objectively evaluate any receiving structure that may be contrived. 

Our aim in this investigation was to determine the ultimate limits above which any 
estimate of the phase or amplitude of the spatial Fourier transform of the object must 
lie. Various receiving structures that approach these limits were found, discussed in 
terms of practical realization, and compared with existing techniques. For our 
purposes, the Fourier transform of the object is more“ convenient to consider than the 


1 



equivalent angular intensity function, since the factors limiting the resolution are more 
easily understood. 

1.2 REVIEW OF THE PROBLEM 

In order to acquaint the reader with the basic resolution problem, we shall review 
the theory of telescope imaging, particularly as resolution is affected by the atmo- 
spheric turbulence (see Fig. 1). 



Fig. 1. Receiver geometry. 

Points on the object are denoted by the angular coordinates _z = (z^, z^), measured 
from the center of the receiving coordinates x = (x^, x 2 >. The light radiated from the 
object is assumed to be noiselike. . That is, similar to black-body radiation, which is 
(i) spatially incoherent across the object, each surface element radiating indepen- 
dently, and (ii) temporally incoherent, with a smoothly varying spectrum. Objects that 
are self- illuminated, as, for example, the sun and other stars, immediately justify this 
assumption. The model is also a reasonable description of most bodies that reflect light 
of natural origin. Indeed, even the illumination of a reasonably rough object by a multi- 
mode laser with bandwidth of a few megahertz can be considered. 

The imaging equations simplify if we assume that the light is almost mono- 
chromatic, and describe the electric field as a narrow-band random process 
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The spatial distribution of the average intensity of this i^adiation constitutes the pat- 
tern that is to be imaged. We assume that the object is stable or can be pe rfectly 
tracked, and describe its intensity by the function 6(z) which is proportional to | E(z,t)|^, 
where the bar is an ensemble average. 

An equivalent description of the object, and one more responsive to resolution anal- 
ysis, is the spatial Fourier transform 

®(k) = y dz 0(z) - (1) 

object 

The vector coordinates k = (k^k^) are the two-dimensional spatial frequencies. By 
good angular resolution in an image, we mean a reconstruction,, in which all spatial fre- 
quency components out to some desired limit appear in the picture within some accept- 
able tolerance. 

1. 2. 1 Image Formation in the Free-Space Environment 

In the absence of any turbulence, the complex envelope of the received field behaves 
as 

E{x,t)oc J E(z,t)e -:)2ir -'-/ K dz, (2) 

obj ect 

where we have assumed (i) that the object is at sufficient distance relative to the aperture 
dimensions that a plane wave approximation is valid, and (ii) that the object subtends a 
small enough angle that sin z « z. 

If an aberration- free lens is placed in the aperture (x plane), the field appearing in 
the focal plane has the form 

E(u, t) °c y E(x,y) w(x) e s/(^) dx, (3) 

where L is the focal length of the lens, and w(x) is a transmittance function describing 
the shading of the lens (see Fig. 2). 

A photograph, like most other detectors, responds to a time- averaged value of inten- 
sity, rather than to the instantaneous value of the incident field. This focal plane inten- 
sity we describe as 

I(U) «<|E(u, t)| 2 > 

cc ^ dxdx' w<x) w(x' ) <E(x, t)]l<x' , t) > 2lT - * ) / XL, J (4) 

where the result of the time averaging is denoted by angular brackets. In general, this 
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does not give a deterministic result because of the random nature of the original radia- 
tion. Practically, however, the bandwidth of even the most monochromatic natural light 
is many times greater than the detector bandwidth. This bandwidth is the reciprocal of 



Fig. 2. Imaging with a lens. 

exposure time for a photograph. Thus the output of the detector represents a sum of a 
large number of independent samples of the radiation process; hence, there is an over- 
whelming tendency for the sample average to approach the expectation of an ensemble of 
source radiation processes. The bracketed term in Eq. 4 is then closely approximated 
by the ensemble average mutual coherence function 

T{x, x 1 ) = E(x, t) !(x\t) « <E(x, t)I§{x', t)). • (5) 

When the object is spatially incoherent, with the aid of (1) and (2) we obtain the 
well-known result 

r<x, x») cc ©[(x-x’)A]. (6) 

That is, the mutual coherence function of the received field is dependent only upon 
the displacement x - x‘> and is simply proportional to the object Fourier transform 
evaluated at k= (x-x 1 )/\. 

The image formed by the telescope can now be obtained. Rather than study the 
pattern in the focal plane, however, it is more convenient to evaluate the image Fourier 
transform 

jT(k) = j* I(-u) e" 3 * 1 *’ 2/ L d u, (7) 
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where the argument of I has a sign reversal to account for the familiar image inversion 
that occurs with a single lens. This image transform becomes simply 


(k) oc 0(k) y dxw(x) w(x+Xk). 


( 8 ) 


1. 2. 2 Diffraction-Limited Modulation Transfer Function 

Equation 8 shows that the spatial frequency components in the object and image are 
directly related by a function 

H(k) = y dx w(x) w(x-fXk), (9) 

which is commonly referred to as the modulation transfer function (MTF). Generally, 
aperture points are either completely transparent or opaque; the pupil function then has 
only two values, zero and one. H(k) is then the area in the available aperture over which 
the mutual coherence function is integrated for a displacement \k. Equivalently, it is 
the common area between the aperture and a replica translated by \k (see Fig. 3). 



The common circular aperture with diameter D has a radially symmetric MTF with 
the general shape indicated in Fig. 4. The upper limit of spatial frequencies appearing 


H(k) 



k 


5 



in the image is then D/X. As an example, a 10- cm aperture operating in the middle of 
the visible spectrum (0. 55 p), has an upper resolution limit of 0. 18 lines/prad. This 
limit is commonly referred to as the diffraction limit of the system, since it is due to 
the fact that diffraction of light from the edges of the aperture diffuses the image. 


I(u) 



In the focal plane itself, this diffraction effect transforms a point source (incident 
plane wave) into an Airy disk surrounded by fainter concentric rings. The general form 
of the intensity is sketched in Fig. 5, and is given precisely as 


I(U) = 


J j [2irD I u | /(XL )] 

L ttD I u | / (XL>) J 


( 10 ) 


Arbitrary resolution criteria have been defined in terms of this point smearing. 

Rayleigh's criterion states that two point sources can be "resolved” when the center of 

the one disk is beyond the first null of the other. This requires that the points be sepa- 

2 , 

rated by more than 0. 61 \/D rad (see Born and Wolf ). 

Although such criteria are reasonable rules of thumb when it comes to visual obser- 
vations, they are rather meaningless in a precise analysis, since they evade entirely 
the effects of background noise in the image. Indeed, if there is absolutely no noise in 
the image, the choice between the hypotheses of one or two sources can be perfectly 
resolved, regardless of the aperture size. In the language of the communication theo- 
rist, it is a singular detection problem. 

It is absolutely necessary, then, in speaking of resolution to include the effects of 
the noise. In any recording device the induced noise is ultimately due to the basic quan- 
tum nature of the photodetection process. With a photograph the noise is associated with 
the individual grains that are developed. In the absence of reciprocity failures and other 
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nonlinear effects, the consequences of the noise can usually be made arbitrarily small 
by increasing the amount of light energy, either by increasing the area of the aperture 
or the exposure time. The diffraction- limited MTF will be maintained, however, only 
when the incoming wavefronts remain undisturbed. 


1. 2.3 Atmospheric Modulation Transfer Function 

The effect of the atmosphere is to reduce the wavefront coherence to only a few cen- 
timeters. The time rate of change of these disturbances is usually measured in milli- 
seconds. These atmospheric effects are commonly described by the introduction of a 
complex multiplicative process, e a ^-’ £,tHjb(x, t)^ eac j l incident wave. The mutual 
coherence function of the received process is now not only dependent upon the displace- 
ment x-x* but also upon the time t and point of origin x. in the aperture. 

Equation 6 for the mutual coherence function now becomes 

r(x,x',t) cc ^ dz 0(z) e _ j 2irs * te-S')A e a(x,z,t)+a(x',z,t)+i[b^,z,t)-b(x , ,z,t)] 

As a function of the displacement x - x ! > r no longer has a one-to-one correspon- 
dence with the spatial frequency components in the object, since the angular variations 
in the disturbances induce a dispersion of spatial frequencies between object and image. 
The analysis is simplified, however, if we suppose that the object is small enough to be 
completely within an isoplanatic patch, that is, a region over which the wave disturban- 
ces are invariant with angle. Equation 11 then reduces to 


T{x, x 1 , t) oc © 


x-x* 
* k 


g a{x, t)+a(x', t)+j[b(x, t)-b(x' , t)] 


( 12 ) 


Although there is now no dispersion of frequencies, the factor of proportionality 
between the mutual coherence function and the object transform is a random process, 
which we denote by 

h<x,x',t) = e a (S» t)+a(x', t)+j[b(x, t)-b(x , ) t)] (13) 

In the focal plane of a telescope, the image is now a time- variant random process. 

In the Fourier transform domain it behaves as 

^ (kit) 01 y w(x) w{x+\k) r(x,x ! ,t) dx oc ©(k) ^ w(x) w(x+kk) h(x, x+kk, t) dx. (14) 

Thus there is a random MTF relating the object and image. When the aperture is 
large or the exposure time long, we expect that the sample MTF will approach the same 
value as an average over an ensemble of atmospheres. Under such an averaging, the 
effects of the aperture dimensions and the atmosphere are separable, and we have 

jT(k, t) oc 0(k) H(k) H A (k), (15) 
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where H A {k) is the average atmospheric MTF. This function has been estimated ana- 

lytically in several ways. Hufnagel and Stanley -5 obtained an approximate solution of the 

wave equation that must be satisfied by the mutual coherence function as averaged over 

4 

an ensemble of atmospheres. Fried has used an approach that gives similar results, 

• 5 

and this is the tack that we follow. Since Tatarski and others have already shown that 
a and b behave as jointly Gaussian random processes, h(x, x',t) can be averaged 
directly to give 


H A (k) = h(x, x+Xk, t) 

-1/2 [D a (|\k|)+D b <|Xk|}] 


- e 


( 16 ) 


Both a and b are assumed locally stationary and homogeneous with structure functions 

D (x-x') = [a(x, t)-a{x’ , t)] 2 (17) 

a 


D b (x-x r ) - [b(x, t)-b(x',t)] . (18) 

The mean value of b is taken as zero but, in order to maintain conservation of energy 



lines/milliradian 

/ 

Fig. 6. Typical atmospheric MTF (Djurle and Back )« 
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regardless of the severity of the fading, the mean and variance of a are related by 


a=-[a-a] 2 . ' (19) 

Because of the exponential behavior of (16), the attenuation is very strong once [ Xk[ 
exceeds the correlation distance. Djurle and Back have made some measurements of 
atmospheric MTF which are representative of the image degradation to be expected. 
These results are reproduced in Fig. 6, and it can be seen that an attenuation of more 
than 20 dB (factor of 10 in H^) is encountered for spatial frequencies greater than 0. 03 
lines/ixrad. Recall that 0. 03 lines/prad is well within the diffraction limit of a 10-cm 
diameter aperture. 

At larger values of displacement the attenuation is even greater, and is difficult to 

measure directly from .a distorted image. In Section IV, some experimental results are 

2 

presented in which it was found that a phase variance of 50 (rad) was common on a long 
horizontal path. Such a variance corresponds to an MTF attenuation of more than 
200 dB. 

1. 2.4 Image Reconstruction by Linear Filtering 

One way in which some of the detail in a photograph can be restored is through 
spatial frequency enhancement techniques. The image from a long-exposure photograph 
is Fourier- transformed, divided by the average MTF, and transformed back, into a 

7 

high-quality picture. Harris used such a technique to make distorted images more 
distinguishable. This technique is similar to that used for edge enhancement in 
restoring telemetered pictures. Its success in the atmospheric problem is rather 
limited, mainly because of the following considerations: 

1. the ever-present background noise is enhanced, together with the signal and may 
tend to dominate the final picture 

2. the enhancement ratio may even be more than the linear dynamic range of the 
original recording 

3. the residual uncertainty between the sample averaged MTF and the ensemble 
value generates an additional error which will also be enhanced. 

o 

Helstrom considered the optimizing of the linear filtering procedure in order to 
minimize the total mean- square error in the picture. This is the optimum filtering 

Q 

technique of Wiener and Lee and requires a priori knowledge of the average spatial 
spectrum of the class of objects considered. 

In many imaging problems an ensemble description is not justified and, indeed, when 
some a priori assumptions are fed to such an 11 optimum 11 processor, the tendency will 
be, especially when the noise is strong, for the estimate to be biased toward the one 
which the processor believes is most likely a priori. For such reasons, it does not 
appear that the high- resolution imaging problem is really too suitable for the Wiener 
technique. 
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1. 2. 5 Image Motion Tracking 
4 10 

Fried and Heidbreder determined the improvements that are realized when the 
random motion of the image is eliminated. Over small values of displacement, the phase 
structure is believed to depend upon the 5/3 power of the separation. With such a struc- 
ture function, it can he shown that over small apertures, at least, the wavefront distor- 
tion is primarily a change in the local angle of arrival. This random wavefront tilting 
appears in the focal plane as a movement of the image. It is interesting to note that a 
5/3 behavior is not far from a 6/3 = 2 behavior, for which the distortion would be 
entirely wavefront tilting. 

If this motion can be eliminated by tracking of the tilt, the resultant average MTF 
can show possible improvement of three times the resolution bandwidth over the uncom- 
pensated case. This resolution bandwidth is defined as the first moment, or center of 
gravity of the MTF. 

This maximum improvement was found to occur when the aperture is stopped down 
to approximately 3 or 4 correlation diameters. For high spatial frequencies, which 
require considerably larger apertures just to overcome the diffraction limitations, the 
reduction in attenuation is unspectacular, image motion being only a small part of the 
total distortion. Moreover, the object must be bright enough to generate a reliable 
tracking signal. 

1.3 OUTLINE OF THE REPORT 

Our intention has been to bring a more fundamental communication theory approach 
to the problem of imaging in a turbulent environment. 

The classical statistical estimation problem is first formulated in Section II where 
we suppose, at least cone optionally, that the incident electric field is the input to the 
signal processor. From this initial formulation we derive the sufficient statistics for 
the incident process. These are the numbers containing all of the information relevant 
to an estimation of the object, and those that a receiver should extract if it is to utilize 
all available information. The necessary measurements are found to be the mutual 
coherence process as appropriately sampled in optical frequency, time, position in the 
aperture, and displacement. Implicit in the sampling operation is a prior local inte- 
gration of the coherence process over corresponding correlation regions of the above- 
mentioned variables. The sample averaging is then shown to give essentially the same 
result as an expectation over an ensemble of source radiation processes, the result of 
which is still dependent upon the instantaneous atmospheric disturbances. The physical 
measurement of mutual coherence then follows, and we explore various interferometer 
configurations of lenses, mirrors, 'etc. terminated with energy detectors. A model for 
the noise that must inevitably accompany the detection is then developed. 

The ultimate performance limitations are the subject of Section III. With the aid of 
the Cramer- Rao inequality we are able to underbound the variance of any unbiased 
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estimate of either the phase or the amplitude of the spatial Fourier transform. Various 
processing structures that approach these bounds are then studied and compared with 
conventional techniques. We find that if the signal- to -noise ratio is very high, optimum 
estimation of the parameters is realized by tracking the phase and amplitude distur- 
bances in the mutual coherence function. At lower values of signal-to-noise ratio, a 
reliable tracking signal cannot be obtained, and a square-law or crosscorrelation pro- 
cessing of the samples appears to be most favorable. 

Finally, with extremely small signal levels it is found that a straight integration of 

the samples may be most advantageous. That is, the techniques of long-exposure 

recording should be employed, perhaps with a device other than a photograph for 

recording, because of effects of reciprocity failure and other nonlinearities. Finally, 

the various proposed processing structures are compared with techniques employed at 

1112 

the present time, including the Hanbury Brown and Twiss 9 intensity interferometry. 

Some results of a parallel experimental investigation are presented in Section IV. 
Photographic and photoelectric measurements illustrate both qualitatively and quantita- 
tively the major effects of the atmosphere as they affect image estimation. 

Derivations that would have required more than a few lines in the main text are rele- 
gated to appendices. 

Only the more recent journal articles are referenced in their original form. Much 
of the older source material has by this time been chewed over by many authors and been 
republished in more digestible forms. For this reason, we have referenced publications 
that are readily available and whose explanations are most suitable for our analysis. The 
reader interested in more source material will usually find that the cited authors have 
provided a plethora of references. 
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II. CHANNEL MODEL 


We have given a rather abbreviated analysis of the general features of the problem 
of imaging through a turbulent atmosphere. This analysis was primarily directed toward 
the image degradation experienced with photographs taken through telescopes. We now 
want to develop a broader model for the problem, and to find the optimal set of field 
measurements that should be recorded and their statistical form. The study of optimum 
schemes of processing these measurements and the performance limitations will be con- 
sidered in Section IV. 

We shall begin with a classical approach to the estimation problem and consider the 
action of a statistically optimum receiver structure when it has complete specification 
of the incident electric field process. By this means, we are able to show that samples 
of the spatial mutual coherence properties of the field constitute a statistically sufficient 
set of measurements. The relationship of these measurements to the average field inten- 
sity pattern across the object will then be discussed. We shall then show that the local 
coherence measurements are related to the spatial Fourier transform of the object 
through a random multiplicative term attributable to the atmosphere and an additive shot- 
noise term introduced by the detector. 

2. 1 RECEIVED ELECTRIC FIELD 

We imagine that the object and our receiving aperture are orientated as in Fig. 7. 

We model the object as an unknown distribution of point radiators, and restrict ourselves 
to that class of objects whose radiation is noiselike. This implies that the radiators are 
spatially and temporally incoherent. By these terms we mean that each element emits 
a statistically independent electric-field process with a smoothly varying spectrum. By 
smoothly varying we mean that significant changes do not occur on a scale of less than 
a few megahertz at least. 

Objects that are self- illuminated immediately satisfy these conditions, since each 
atom is excited independently and decays with a time constant that is usually very short. 
Only when the radiators are coupled, as in a resonant laser cavity, must we be concerned 
with the temporal characteristics of the radiation. 

When the object is a secondary scatterer, we must be somewhat more careful 
because if the illuminating light is either highly temporally or spatially coherent, the 
light received from the object may still carry some of these properties. Often, however, 
when the surface of the object is rough and has some vibrational motion, these coherence 
properties are destroyed. What we must exclude from our model are such phenomena 
as the 11 grainy” pattern sometimes observed on an object illuminated by a laser beam. 

The function that is of interest to us and constitutes the picture of the object is the 
average intensity profile which we denote 0 (z), where the subscript v allows for a pos- 
sible frequency dependence of the radiation. An equivalent description of the object, 
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i>2 Object 



Fig, 7. Problem geometry. 


and one that is often more responsive to resolution analysis is the spatial Fourier trans- 
form, which we denote 0 (k), where the vector k represents the two-dimensional 
spatial frequencies. 

The electric field radiated from the object is modelled as a Gaussian (normal) ran- 
dom process. This process is most conveniently described with a Fourier series rep- 
resentation. We first subdivide the entire observation time interval T into segments 

of length AT, where AT is the correlation time of the atmospheric disturbances. That 

th 

is, we can consider the disturbances as constant over such a period. Over the 1 such 
interval the field emitted from point p on the object is expressed as a sum of frequency 
components as 



( 20 ) 


( 21 ) 
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Only those frequencies v = m/AT that fall within the optical spectral range of the 
object give significant contributions to Eq. 20. We assume that the field strength E^(t) 
is referenced to unit distance from the element. In the absence of an inhomogeneous 
atmosphere; the corresponding complex representation of the field received over the 
aperture coordinates x = (x^x^) * s 




P’ 1 ’— -jZTTj-x/tXR) 
R 


( 22 ) 


where | = (£ , £ ) are the coordinates of the point p, and X is the wavelength associated 
with frequency v. 

Two simplifications are implicit in arriving at (22). First, we have assumed that the 
aperture is sufficiently small compared with the radius R of the incident spherical wave, 
so that a plane wave approximation is valid. Second, we have ignored the average prop- 
agation delay, since only the relative delay over the aperture is of importance. 

5 

When turbulence is present, we follow the lead of Tatar ski and others and introduce 
a random amplitude and phase term of the form 


a (x, i, t)+jb (x, £, t) 

e v ' v . (23) 

Theoretical and experimental results strongly suggest that both a and b are well 
modelled as Gaussian random processes. The results here do not require this Gaussian 
hypothesis, and only the correlation scales are needed. The correlation times typically 
have been found to be of the order of 1 msec or more with a spatial scale of, perhaps, 
a few centimeters. The strength of the fading is described by the variance of a, which 
has been found to be as much* as one at times. The variance of the phase process is often 
extremely large, but in our work the absolute phase is unimportant. Rather the phase 
difference between aperture points is of more significance, and has been observed to 
vary many radians. This spatial phase incoherence is the most troublesome as far as 
imaging is concerned, but, as we shall see, the amplitude disturbances also play an 
important role in the performance of estimator structures. 

The description of the received field is not complete without including some back- 
ground noise. This may be sunlight, moonlight, and so forth scattered from particles 
in the air. We simply represent it as an additional Gaussian component n^ m (x) in the 
received field. 

With the effects of the atmosphere and the noise included, (22) now becomes 


, , V Pp ’ l ’ m H2ir£-x/(XR) V2S.l*t)+jb v (x4,t) 

^1, m-^ R e "" e +n l,m<2>- 


(24) 
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2.2 SUFFICIENT FIELD MEASUREMENTS 


The functions ^(x), though a complete representation of the received random pro- 
cess, cannot be practically measured if for no other reason than that it would require 
so many functions to represent the process over any reasonable interval. We now seek 
a simpler set of measurements that still retain all of the information relevant to the 
intensity pattern on the object. We continue to treat the problem classically and ask 
what a statistically optimum processor would do with the functions E, (x) if it actually 

1 y xn 

had them available. Such a receiver extracts all of the relevant information, by calcu- 
lating a measure of the probability of the received process conditioned on all acceptable 
hypotheses. This conditional probability is very difficult to express directly, since the 
processes E^ (x) are products of normal and log-normal forms. We circumvent this 
difficulty by conditioning the probability function on the atmospheric processes a and b, 
as well as on the unknown object pattern 0. The processes m (x) are then Gaussian 
under the additional condition, and the probability function is well known. The condi- 
tioning on the atmosphere can be removed by averaging the result over a and b; however, 
this last operation is unnecessary in the determination of the simpler yet sufficient mea- 
surements. _ 

As we have mentioned, the spectra of our objects are smoothly varying on a scale of 
at least many megahertz. Thus when we choose AT as 1 msec or so we expect the set 
of spatial random processes m (x) to be independent in the disjoint time intervals 1 
and orthogonal frequencies m. The probability of the entire set of processes is then 
simply a product of probabilities of each term taken separately. 

Before we can express the probability of a particular sample function we must reduce 
to a random -variable representation. Such a representation is most conveniently 
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obtained with a Karhunen-Loeve expansion (see Van Trees ). We first require the spa- 
tial correlation properties of the received process. These properties are completely 
contained in the two complex correlation functions 


1, m — — 1, m 


(x) fe. (x’) 

— 1, m — 


(25) 





(x,x’)=E. (x) E. (x'). 

m - - 1, m - 1, m — ' 


(26) 


Under the natural assumption that the complex coefficients describing the object and 
noise radiation have uniformly distributed phases, (26) is identically zero. Also, the 
real and imaginary components of each sample function turn out to be uncorrelated. Each 
of these components then has the same correlation function 

R l, m ( — ’ — } = Re [E l,m ( ^ Re [Ei.mte’)! 3 Im E B i f m<2>] Im [E l,m ( ^ )] 

= 2 Re t r i.m (s ’ £,) l {27) 
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The processes E^ m (x) are now expanded in the orthogonal series 
00 

E, (x) = ) E, ij>, (x), 

1, m'-' X, m, n T l, m,n v - / ’ 


n=l 


(28) 


where the ^(x) are real functions that have been chosen as the orthonormal eigen- 

functions of the equation 


c, ik • (x) = 

lj nij n I) m, n 


\ ^(x') R, (x,x’) dx’. 

J L, m, n — 1, m — ' — 

aperture 


The coefficients in the expansion (29) are determined by 


(29) 


E 


1, m, n 


I 


E, (x) v| j (x) dx. 

1, m - 1, m, n - 


aperture 


(30) 


These coefficients have independent zero-mean real and imaginary components, each 
with a variance c^ , the corresponding eigenvalue of (29). 

If we approximate (28) with the first N terms, the conditional probability then is 


P<E lf m (x)/ e > a,b) 



( 31 ) 


Multiplying the results for all time intervals and frequencies, taking the logarithm, 
and retaining only the part that is dependent upon the electric field leaves 


T/AT 


00 


III 


(32) 


where the summation on n has now been extended to +». This triple summation (32) is 
then the operation that an optimal receiver would perform on the electric field, but the 
receiver is required to know the c^ for each possible hypothesis concerning 0, a, 

and b. 

The determination of a sufficient set of measurements comes, as we now show that 
the summations in (32) can be reduced considerably even with 0, a, and b unknown. 

We first rewrite (32) as 


T/AT w 

I 2 E l. m® K *’> ' 


1=0 m=0 


(33) 
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where 


«i. »<*•*> = i 


00 (x) iLl (x f ) 

— 1, m, n v - ; T l, m, n - J 


2c 


n= 1 


1, m, n 


(34) 


is the inverse kernel of the correlation function R. (x, x 1 ). 

i f m 

Now we can safely assume that over some range of the frequency index m, the kernel 
Q will remain virtually unchanged. We then rewrite the summation on m in (33) as 


M(q+1)— 1 

I < > = 1 Q l,M„te2'> I E l,m® 


(35) 


m=0 q=0 m=Mq 

and with the aid of Parseval's theorem we find that expression (33) is proportional to 

(1+1) AT 


T/AT M 
| di dx' ^ £ %, M qteS') 


E (x, t) (x 1 , t) dt. 


(36) 


1=0 q=0 


l'AT 


The term E v (x,t) in expression (36) is just the narrow -band field process resulting from 
the M Fourier frequency components, and can be obtained by bandpass -filtering the inci- 
dent field between frequencies Mq/AT and M(q+1)/AT, v « Mq/AT. 

The integration in (36) is then simply a sample time average of the coherence between 

the fields at points x and x 1 within the optical band about frequency v. This average is 

\ 

commonly referred to as the mutual ^ coherence process. We use the word process rather 
than function to denote that it is still a local sample average, rather than an ensemble 
average. 

The number of terms M that can be brought within the partial sum of (35) is an 
important parameter of the measurement, as it sets the optical bandwidths of the set of 
filters that the receiver should employ. This bandwidth we denote 

Av = M/AT. (37) 


We can estimate the width of this band by considering the correlation function 


r i, m ^x')=E ltm {x)ii >m (x-) 


= 1 


l p P,l, 


m' 


• (x-x f )/(VR) 
e 


h (x, x 1 , £,t) + n, (x) n, (x r ), 
v “ — — 1 1, m - 1, m'— ; 


( 38 ) 


where h is the atmospheric function defined as 
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( 39 ) 


a v (2£> 1- t)+a (x 1 , i, t}+j[b (x, i, t)-b (x 1 , i, t) ] 
h v (x,x',|,t) = e v v " v v “ - 


The dependence upon m in the terms |p , 1 2 and n, (x) n, (x 1 ) describes the 

1 p, 1, m ' 1, m- 1, m- 

spectral behavior of the object and noise background. We shall suppose that these are 
sufficiently wideband that the limit on Av is due to the behavior of the terms 

(x-x')/(\R) 

e (40) 


and 


h v (2»2£'’ *)• (41) 

To see the general order of Av to be encountered, consider the following example: 

1. maximum field of view 

|||max/R = 100 prad (~25 seconds of arc) 

2. maximum aperture displacement 
|x-x T I " 1 meter. 

We find that the phase angle associated with (40) changes only about 1 rad when the 
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frequency is swept over some 0. 5 X 10 Hz. In terms of wavelength bandwidth this is 

o 

almost 50 A for light near the middle of the visible spectrum. 

The optical correlation bandwidth of the term (41) is also significant in determining 
the allowable Av. Very little experimental work has been reported on the frequency 
dependence of a and b, but we can still get a rough estimate with a ray theory argument. 
Under the assumption that the index of refraction of the air itself is constant over the 
small optical bands of interest, the ray analysis predicts an attenuation term a inde- 
pendent of frequency, and a phase that is linearly related to frequency. Thus if at one 
particular frequency the phase difference between aperture points is 10 rad, a frequency 
10% higher would show a corresponding phase difference of 11 rad. Admittedly, a ray 
analysis is weak on a long path, diffraction effects being important. Nevertheless, we 
can still convince ourselves that the optical correlation bandwidth Av is of the same gen- 
eral order as that predicted in the previous paragraph. 

Although interpretation of expression (36) gives the general form of the sufficient 
measurements, namely time and frequency sampled coherence measurements, we still 
seek some additional reduction. Rather than have to measure this coherence process 
for all continuous pairs of points x and x‘> we want to get a spatial sampled measure- 
ment. Showing that this reduction is possible requires demonstrating that the kernel Q 
has characteristic scale dimensions x and x. 1 over which there is little significant change. 

This is most easily demonstrated if we look at a special case. Let us suppose that 
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the background noise power is high compared with the signal, and that the noise corre- 
lation function is impulsive in the displacement variable x*~x f , or at least has a correla- 
tion scale shorter than that of the field from the object. As we shall see, this essentially 
means that the noise power appears evenly distributed over a region of the sky larger 

than that subtended by the object. Under these conditions, we can invoke the CUP results 
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(Coherently Unestimable Parameters — see Van Trees ). We are able to approximate 
the kernel 


m m 


where we have supposed that for the stationary noise process 


(42) 


n, (x) n, (x 1 ) = N. (i (x-x'). 
1, m 1, m - m o ’ 


(43) 


Expression (17) then reduces to 


t/at ^ 

I l 

1=0 q=0 


N. 


(1+1) AT 

“ \ dx \ |E (x,t)| 2 dt 

V/Tn O J V 

t=lAT 


Mq 


N 


f- I dX I dx ’ E l,Mq feS) 


(1+1) AT 


Mq 


r 


E v (x,t) E v (x’,t) dt 


t=IAT 


(44) 


The first term in (44) is simply a measure of the total incident flux in each band, 
and is an almost trivial measurement. In the second term the kernel is now the corre- 
lation function whose dependence upon x and x 1 is known in terms of the atmospheric 
function h. Consider then just the integral part of the ^second term of (44). Changing 
variables so that x 1 = x + y, we can rewrite it 

(1+1) AT 

J dy J dx R 1Mq (x, x+y) j E v (x,t) & v (x+y.t) dt. (45) 

t=lAT 


Suppose now that we divide up the aperture into cells, each of side Ax, where Ax is 
the x. correlation scale of the atmospheric function h. Numbering these cells with the 
index r = 1, 2 . . . , we can then rewrite (45) as 


th .. 
r cell 


(1+1) AT 

f d Z £ R l,Mq ( W£> I d ^ J VS’ t)fi v ( SV’ t)dt 

t=lAT 


(46) 
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Thus we find that sufficient measurements can be obtained by a local integration of the 
mutual coherence in both time and aperture space. 

We still have the displacement variable y which can also be quantized. We can see 
from Eqs. 36 and 39 that the dependence upon y is in the atmospheric function h and 
the term 

~j2v£- (x-x')/(XR) +j2*rry £/(XR) 

e = e “ " . (47) 

The correlation scale for y in the atmospheric function is the same as for the x vari- 
able. It is usually the second function (47) which sets the quantization scale of y for 
reducing (46) further, and it is the field of view subtended by the object that is the deter- 
mining factor. We require that the exponent term 2iry * £/(XR) change less than approx- 
imately 1 rad over the local integration on y. Indeed, the operation of integrating the 
displacement variable over a Ay X Ay region is really equivalent to restricting the field 
of view of the receiver to a solid angle of 

S2 = \ 2 /Ay 2 . (48) 

This can be thought of as spatial filtering to suppress background noise. 

In summary, a statistically sufficient set of data can be obtained by local measure- 
ment of the mutual coherence process denoted 

r v (x, x+y, t) = < E v <x, t) ^ v (x+y, t) ), (49) 

where the angular brackets represent the local integrations in 

1. time over the correlation time AT 

2. frequency over the correlated optical bandwidth Av 

2 

3. position over the correlation region Ax 

2 2 

4. displacement over the region Ay = X /ft. 

A further point should be mentioned before passing to the next topic. Throughout 
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the foregoing analysis we treated the electric field as a scalar quantity. Saleh has 
shown that the atmosphere introduces no significant depolarization effects. If the inten- 
sity of the emitted radiation is the same for both polarization components, and these 
components are uncorrelated (uniform polarization angle), the mutual- coherence samples 
for each component should just be added together to form the sufficient measurements. 

2, 3 BEHAVIOR OF SAMPLED MUTUAL-COHERENCE PROCESS 

We now want to consider the statistical properties of the measurements defined by 
Eq. 30. There are basically three random phenomena reflected in the data: 

1. the radiation processes (both object and background) 

2. the atmospheric turbulence 

3. the noise disturbances introduced by the detector. 
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The first of these phenomena is most easy to treat in the optical case and, in fact, 
is almost negligible compared with the other two. This is because the large number of 
essentially independent samples of. the process entering into each sample integration 
lets us substitute the ensemble average as the result. 

2. 3. 1 Ensemble -Average Mutual -Coherence Process 

If we average over the entire ensemble of object and background radiation processes, 
the average data samples are proportional to the correlation result of (38), where we 
still treat the atmospheric disturbances as parameters rather than random processes, 

r v (x, x+y, t) cc r i m (x, x+y) 

In I 2 

^P> 1, m -j2ir£-y/(\R) * 

2 e h v ( -’ + n l, m W n l,m ( - + l)’ ( 50 ) 

P R 

where 

£ is the location of object point p 
1 is the time index 
m is the frequency index 
v is the frequency = m/AT 
\ is the wavelength = velocity of light/v 
h is the atmospheric disturbance (Eq. 20) 
n l m iS bac kground noise process 

jpp ^ m is proportional to the average energy radiated from the object (for sim- 
plicity, we suppress the factor 6^/4 which really relates the energy to the mean- 
square value of field coefficients). 

Representing the object as a spatially continuous collection of radiators allows us 
to replace the summation in (50) by an integral 




P 


r d| ~ 

> - J rt Vi )( >• 


( 51 ) 


^ th 

where 0 (£) is the average energy flux received on the m frequency during a AT time 

interval from an area d| = d^d^ on the object. Time dependence (subscript 1), has 

been dropped in 0, since we take for granted that the pattern on the object is stationary. 


We now replace the linear coordinates by angular variables 



( 52 ) 
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and 


is 

Z 2 " R ' 


( 53 ) 


It is assumed that the object is of small enough angular size that sin z ~ z. 
Expression (50) can now be written 

7 n +jZnz-y/X “ 

r v (x,x+y,t) °c J dz 6 v (z) e \(x,x+y, z,t) + n^ m (x) n l m (x+y). (54) 

6> v (z) is now proportional to the energy per solid angle and has the units: watts/(unit 
bandwidth), (unit solid angle), (unit area on aperture). 

If we consider, for a moment, the special case of a free-space environment (h = 1), 
expression (54) reduces to simply 


r v (x.x+y,t) cc © v (-yA) + n^ m (x) n^fx+y), (55) 

where © is the spatial Fourier transform of the object defined as 

J -jZttz * k 

dz 6 v (z) e (56) 

We can look at the background -noise process as being due to a continuous distribution 
of radiators on a shell centered on the aperture. If we denote this distribution by the 
function N^fz), assuming that the received noise is both stationary and homogeneous, 
the spatial correlation then follows from analogy with (55) and (56) as 

; 7 r +j2trz • y/X 

\ m ( *> n l, = J N v (z) e < 57 > 

Because of the filtering operation in the receiver structure (local integration on y 
restricting the field of view}, we need only be concerned with the nature of N (z) in the 
neighborhood of the object. Usually the noise will appear uniform over a reasonably 
large area, and we can simply denote the spectral flux density of the noise radiation as 
N v watts/ (unit bandwidth), (unit solid angle), (unit area on aperture). 

The ensemble average value of our measurements, which contain a local integration 
on t, v, x, and y becomes 

7 2 2 C +j2irz.y/\ „ 

r v (x.x+y,t) cc AT Av Ax Ay J dzS^zJe " h^x, x+y, |, t) + AT Av Ax N v 6 0 . 

(58) 
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2. 3. 2 Convergence of Sampled Mutual- Coherence Process 
to Ensemble-Average Value 

The foregoing analysis has all been concerned with the ensemble average value of 
the mutual- coherence process. We now want to show that, with high probability, the 
sample average that we obtain by integrating over the channel coherence time AT yields 
essentially the same result as the ensemble average. By ensemble average we mean an 
average over only the random fluctuations in the emission process; we still have depen- 
dence upon the phase and amplitude disturbances introduced by the atmosphere. 

The convergence of the sample average is demonstrated by the mean- square error 

2 

e defined as 



f v (x,x , ,t)-r v (x f x , ,t)| . 


( 62 ) 


This turns out to be (see Appendix A) 

2 4 

C - AT Av Ax 1 


2a (x, t) 

2a (x’.t) 

© v (0) e v 4 N v a 

© v (0) e 4 N v £2 


(63) 


If we normalize this error by dividing by the mean- square value of the mutual coherence 
measurement, we obtain 


2a (x, t) 

® (0) e v + N n 

V V 


2a (x’,t) 

© (0) e v +N n 

V ' V 


Iryx^t) 


2 ” AT Av 


« 2[a (x, t)+a (xSt)] 

l® T [(S-=?)A ]| 2 e r ~ v 


(64) 


This normalized error drops inversely as the product AT Av, the number of essen- 
tially independent samples of the emission process enclosed in the sample average. The 
other factor is often defined as a "fringe visibility function" (FVF). 


FVF = 


. I? ^[a (x,t)+a (x’,t)] 


2a (x, t) 

® (O) e v *' + N n 

v v 


2a (x',t) 

© (O) e v 4' N U 

v v 


(65) 


As an example of the strong convergence of the sample mutual coherence toward the 
ensemble -average value, consider the case when the FVF is, say, 10 . Were it not 

for the noise term in (60), this would mean that the spatial frequency component in which 
we are interested has a squared magnitude only 1% of that of the zero-order component. 

O 

With a conservative time -bandwidth product of 10 , the resultant mean-square error in 

/ 

the sample coherence average is only one part in 10 . The error tends to increase with 
the depth of fading; however, only rarely would the error be more than one or two orders 
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The delta function 6,-. is used to show that the noise contributes to the correlation only 
0,y J 

at zero displacement. For all other discrete sample points of y the noise field gives 
no correlation. 

In evaluating the y integration on the background noise term (38) we assumed that 

the noise power was uniform over a. region considerably larger than that subtended by 

2 2 2 

the object. If we divide expression (58) by Ay , replacing X /Ay by the field of view 
£2, we get the somewhat more agreeable form 


r v (x,x+y,t) cc 


AT Av Ax 


+j2irz • y/X 


J dz BJz) e h Tj .(x, x+y, z, t) 


+ N fiS„ 
v O, y 


(59) 


When there is no atmospheric turbulence, the first term in Eq. 59 reduces simply 
to the Fourier transform of the object evaluated at spatial frequency (k) = -y/x. In gen- 
eral, when turbulence is present, the angular variations in the atmospheric function h 
induce a dispersion of spatial frequencies; that is, there is no longer a one-to-one 
correspondence between spatial frequencies on the object and mutual coherence mea- 
surements on the receiving aperture. We simplify our analysis considerably, how- 
ever, by restricting ourselves to those cases in which the angular dimensions of 
the object are less than the correlation angle of the disturbances. In section 2. 3. 3 
we shall consider the possible applicability of the isoplanatic results to the more 
general case. 

In this isoplanatic case the function h is independent of z, and we then have 


r (x, x+y,t) cc Av AT Ax 2 © v (~y) h^fx, x+y, t) + Av AT Ax 2 N 


6 — , 

v O, y 


( 60 ) 


where 


a (x, t)+a (x+y, t)+j[b (x, t) -b (x+y, t)] 

h^x.x+y.t) *e v v ~ v ' v " . (61) 

We again have a one-to-one correspondence between the Fourier transform © 
and the average mutual coherence, but the multiplying factor h is a random process 
that prevents us from obtaining © easily from single measurements of T\ 

The effect of the background noise, since we have assumed that it has uniform inten- 

A 

sity across the field of view, is simply to raise the zero displacement value of T above 
the @{0) value. We can think of the background radiation as simply an intensity bias 
to the object pattern. 
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of magnitude greater. Because the' error is so small, we can neglect it in comparison 
with the multiplicative disturbances introduced by the atmosphere and the inevitable mea- 
surement noise resulting from the quantum nature of the detectors. This is very impor- 
tant, as it means that not only is the mutual- coherence process a statistically sufficient 
description of the received electromagnetic phenomenon but the original uncertainty in 
the emission process is insignificant in the sampled values. 

For this reason, we shall now discard the circumflex and denote the sample mea- 
surements r v (x,x+y, t), since they are essentially identical to the ensemble values. 

2. 3. 3 Extension to the Nonisoplanatic Case 

Much of the preceding analysis has been predicated upon the condition that the object 
was entirely within an isoplanatic region. That is, both the phase and amplitude distur- 
bances induced by the turbulence were invariant with angle over the dimensions of the 
object. If an object is small enough, this condition will surely hold. Often, however, 
there will be a variation in the angular dependence of the disturbances as we sweep 
across the source and we shall no longer get the simple relationship (Eq. 60) between 
the mutual- coherence process and the object intensity pattern. 

In this more general case we have the result of Eq. 59 (disregarding the noise term) 

p j2vz, y/X 

r v (x,x+y^) = j dz © v (z) e ~ h y (x, x+y> z, t), (66) 

or rewritten in terms of ®(k) 

r C -j2irz -(k-y/X) 

t) = J dk © v (k) J dze h y (x, x+y, z, t), (67) 

where the inner integral determines the severity of the spatial frequency dispersion, 
because of the angular variations of the disturbances. 

We can often reduce this more complicated form to something quite similar to the 
isoplanatic case. Suppose that we divide the field of view into a grid of areas, each of 
which is within an isoplanatic region. Furthermore, suppose that at the receiver we are 
able to separate the mutual coherence process into the separate contributions from each 
of these cells. If we can do this, we can then decompose the channel into a number of 
parallel channels, each of which behaves in the isoplanatic manner, although the partic- 
ular atmospheric disturbance may be different on each channel. We can look at each 
channel independently, and build up a composite picture of the object. 

The problem in such an approach is distinguishing the contributions to the mutual- 
coherence process from each individual. isoplanatic region. Whether we can or cannot 
make this separation depends upon the angular scale of the disturbances relative 
to the correlation scale as a function of the aperture coordinate; that is, whether 
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the aperture coherence area Ax 2 is capable of "resolving" the angular correlation scale. 

The separation of the mutual coherence process into angular components is accom- 
plished by first mathematically Pourier-transforming on the variable y to generate a 
pseudo image 

r a -j2iru • y/\ 

I(u,t) = J dy r v (x,x+y,t) e “ . (68) 


This pseudo image differs from the regular image formed by a lens (see sec- 
tion 1.2.1), in that there is no x integration of T beyond a coherence area. If we ignore 
the background radiation, which only adds a uniform intensity bias, (68) becomes 

a p -j2ir(u-z) -yA 

I(u, t) = \ dz 0 v (z) J dye h v (x, x+y. z> t) (69) 

which can be put in the simpler form 


I(u. t) 


I 


dz e^(z) H v (x,u-z, z, t), 


(70) 


where 


w, z,t) = § dy 


-j 2 ttw • y/\ 


h v (x, x+y, z, t). 


(71) 


It is in this pseudo -image plane u that we try to separate the various isoplanatic 
regions. The effect of the impulse response function H in (70) is to give a spreading 
or smearing of the function 0. This impulse response is a random process that depends 
upon frequency, time, and position coordinate x, but more important, at present, it 
may also depend upon the angular coordinate z. If the change in H with respect to the 
variable z is on a scale that is smaller than the image spreading described by the inte- 
gral operation (70), any point u will be comprised of contributions from the object that 
have suffered different distortions. Thus, to be able to reduce the channel to the desired 
form, we must guarantee that the width of the impulse response H (in the coordinate w 
of Eq. 71) is less than the characteristic angular scale size or angular correlation 
distance Az. The width of H can be estimated from the reciprocal width of its Fourier 
transform, h. We have previously stated that the phase and amplitude disturbances 
remain correlated over a distance Ax. Thus h will show no significant changes with 
y on a scale smaller than Ax, and the width of the impulse response is then roughly 
\/Ax. The pseudo image I(u) can then be broken into isoplanatic regions, over each of 
which the impulse response is invariant with z, provided that 


\ 

Ax Az 


< 1 . 


(72) 
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In radio- scatter communication there is an analogous condition in which the channel 
is said to be underspread, provided that the product of the frequency-dispersion band- 
width and the time spreading is less than one. We recognize that l/Az is the spatial fre- 
quency dispersion, and X/Ax is the angular spreading, so we appropriate the term under 
spread to describe the condition (72), We refer to the converse situation as the over- 
spread channel. 

In the rest of this report we shall consider only the under spread condition. 

2.4 MEASUREMENT OF THE MUTUAL- COHERENCE PROCESS 

We have shown that a sufficient set of statistics can be extracted from the received 
field through sample measurements of the mutual-coherence process 

r v (x, x+y, t) = <E v (x, t) l v (x+y, t) ), (73) 

in which the angular brackets represent a local integration over the optical correlation 
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bandwidth Av, the correlation time AT, and the correlation area Ax within the field of 
view, ft. These sampling parameters are dependent primarily upon corresponding cor- 
relation properties of the turbulence -induced amplitude and phase disturbances. 

We now turn to the physical measurement of these sufficient statistics. At low fre- 
quencies (radio frequencies) the coherence of an electromagnetic field is determined by 
analogue correlation of the voltages induced on two antennas. At optical frequencies, 
however, such direct measurements are precluded, and instead the phenomena of wave 
interference are used. The magnitude and the phase of the mutual -coherence process 
are determined from the degree of constructive or destructive interference obtained when 
the fields at points x and x+y are added. 

Just about all interferometers exploit this basic principle, and differ only in the way 
in which the fields are added together; that is, in the vax'ious forms of temporal and spa- 
tial modulation which are introduced so that the interference effect can be distinguished. 
We shall describe several of the interferometer types that seem suitable for sampling 
the coherence process in the required manner. We assume that the incident light has 
been already filtered to the correlation bandwidth Av, 

2. 4. 1 Simple Interferometer 

Let us consider the simplest scheme for measuring field coherence. From such an 
understanding we can then better appreciate more complicated systems. Probably the 
most basic configuration is the one sketched in -Fig. 8, in which the fields from the two 
arms of the interferometer are superimposed by means of the mirrors and beam 
splitter. 

One of the mirrors is usually tilted slightly by a small angle so that a spa- 
tially varying phase, 2 ttu‘ y/\, exists between the two superimposed field patterns. 
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u plane 


v plane 


The intensity of the total light in either the u or v plane has the form 


I(u,t) 


E v (p=u,t) + E (x'=u+y) e 


+j2iru # v/X 


- ^ r v (u, u, t) +“ r v (u4-y, u*fy, t) + |r v (u,u+y,t) | cos [2im *y/\+<})(u,u+y, t)], (74) 


where cj>(u, u+y, t) is the phase angle of the complex mutual coherence, and factors of 
l/Z are due to the energy division by the beam splitter. By measuring the amplitude 
and spatial phase of the fringes, we obtain the corresponding mutual -coherence mea- 
surement. 

The usual way of extracting the amplitude and phase of the fringes is to correlate 
the spatial pattern with both sine and cosine patterns of the same spatial frequency. A 
direct electrical analogue signal can be obtained by scanning a mosaic of photoelectric 
detectors. This can be an array of separate elements so connected as to give the appro- 
priate correlations, or it can be an internally scanned detector such as an image orthicon 
or vidicon tube. This scanning can be done later if a fast exposure photograph is 
taken of the fringes (exposure time less than the correlation time AT). 

There are also several inexpensive electromechanical means of determining the spa- 
tial phase and amplitude of the fringes. One can simply scan the fringes with a moving 
pinhole in front of the detector surface, and thereby generate a time -variant electrical 
output from which the phase and amplitude can be* determined. This is an inefficient 
means, however, since most of the available light energy is wasted. A more efficient 
scanning mechanism incorporates a glass plate having alternate parallel opaque and 
transparent rulings with spacing matched to that of the fringes, as shown schemat- 
ically in Fig. 9. 
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Fig. 9- Fringe scanner. 


When the plate is moved uniformly in a direction perpendicular to the rulings and the 
fringes, the transmitted light gives a periodic detector output signal. The phase and 
amplitude of this electrical signal are the same as those of the fringes. 

The variation in transmittance on the plate need not be exactly sinusoidal, indeed a 
square wave is fine. The higher spatial harmonics in the ruling pattern only correlate 
with higher spatial frequency components in the fringe pattern. If such components 
exist in the fringe pattern, their effect in the detector output can be suppressed by 
electrical filtering, since they appear as corresponding high-frequency terms in the 
output signal. 

Ramsey and Kobler*^ have cleverly employed this fringe -scanning principle by using 
a rotating disk with radially ruled lines. When the fringes are confined to a small area, 
the rulings appear essentially parallel. A somewhat similar device was employed in our 
measurements (see Section IV). An inexpensive Ronchi ruling {a glass plate with parallel 
opaque and transparent lines) was moved back and forth in reciprocating fashion to scan 
the fringes. 

2 . 4. 2 Fresnel Biprism 

Another way of implementing the same basic interferometer of Fig. 8 (but avoiding 
the mirrors and beam splitters) is by means of a Fresnel biprism whose operation is 
indicated in Fig. 10. 

The light incident on the top half of the prism is deflected downward by an angle y, 
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and that incident on the bottom half deflected upward by the same angle. Over a plane 
L units behind the prism, the intensity is 

I(u, t) = T v (x, x, t) + r v (x\x T ,t) + 2 |r v (x, x 1 , t) | cos (2 ttu 2 sin y/X + 4>(x, x 1 , t)), (75) 

where 

x = u x ,u 2 + L tan y ( 76 ) 

and 

x 1 = u i * u 2 “ L tan Y ' ( 77 ) 

When the angle y is small we can approximate both tan y and sin y by y. 

Fringes appear in the u^ direction, and these fringes have the same phase and ampli- 
tude as those of the mutual -coherence process. The displacement in the measurement 
is related to the distance L. between the prism and the recording, and the effective dis- 
placement is 

y = (0, 21, tan y/X). (78) 

Moving farther back from the prism yields a larger value of spatial frequency, and we 
have the advantage that the fringe spacing remains the same, thereby simplifying the 
scanning device. By also rotating the prism, all values of y can be considered. 

As we have mentioned, the mutual- coherence process should be sampled in a number 
of dimensions. Presumably, the correlated optical frequency bands can be isolated with 
a set of bandpass interference filters, and each detector output can be sampled appro- 
priately in time. Isolation of the measurements to the aperture correlation areas can 
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be realized by locally integrating the fringe patterns, thereby sampling with a different 
detector element for each correlation area. 

2, 4* 3 Lloyd 1 s Mirror 

The Fresnel biprism samples the coherence function in x, but without movement of 
the prism yields only a single value of y* A device which can sample the coherence 
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Fig* 11* Lloyd ! s mirror* 

process in the x^ aperture dimension and the y^ displacement dimension is the Lloyds 
mirror arrangement shown in Fig. 11. The mirror is placed at an angle y to the 
incoming radiation* The intensity of the field on the u plane then is 


I v (u,t) = |E v (x 1 =u 1 ,x 2 =u 2 ,t) +E(x 1 =u 1 ,x 2 =u 2 -2u 2 cos\,t) exp(j2ir(2u 2 siny)A) [ 2 . (79) 


The local interference fringes at (u^, u 2 ) then give us the mutual- coherence measure- 
ment 

r[x-u j, u 2> x+y = Uj,u 2 ( 1-2 cos y)3- ( 80 ) 


That is, the displacement is 

y = (0, -2u 2 cos Y). (81) 

If additional mirrors and prisms are used, the intensity on the measurement plane u 
can be made to take the form 


I v ( u,t) = 


E 


.(x=u,t) + E v (x=-u,t) e 


j2-ry -u/k 


(82) 
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and the fringe producing term is 


2 Ir v (x=u, x+y=~u, t) | cos [Zttv* u/X +<i>(x=u > x+y=-u, t)]. ' (83) 

The local phase and amplitude of this term yield the mutual coherence for all values of 
displacement as ,y = ~2u/X. It should be noted that each value y is sampled with a sep- 
arate pair of aperture points. 

All of these interferometer configurations can be modified by employing a time jitter 
in one of the paths; for example, by twisting or moving a mirror. The phase and ampli- 
tude of the field coherence can then be obtained from the temporal phase and amplitude 
of each detector output, rather than by spatially scanning fringes. 



Fig. 12. Aperture reduction. 

One might not necessarily see these various interferometers in just the form 
described. By various clever uses of special prisms and mirrors, more practically 
useful devices can be constructed that still exploit the same basic ideas. 

It should be pointed out that the prisms and mirrors used in these schemes need not 
be as large as the actual receiving aperture because the field size can be conveniently 
reduced by means of an objective lens (or mirror) and a lens of shorter focal length 
(see Fig. 12). The field scale is then reduced by the ratio of the focal lengths of the 
elements. 

2. 4. 4 The Lens as an Interferometer 

We recall that the mutual- coherence process can be determined by looking at the 
Fourier transform of the image formed by a simple lens (see Fig. 13). 
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2 . 2 
each of area Ax and separated by some distance y (see Pig. 14). With Ax matched 

to the correlation area of the atmospheric disturbances, (86) can then be rewritten 


«^ v (k,t) = Av r v (x o ,x o -Xk,t) ^ w(x) w(x-Xk) dx. 


(87) 


The function defined by the integral in expression (87) determines the values of dis- 
placement and corresponding spatial frequencies which can be extracted from the image 



transform We denote this integral as a transfer function H^-(k) , which is determined 
in exactly the same way as the ordinary diffraction-limited MTF (see section 1. 2. 2). 

H M { - } = J W( - } dx. (88) 

Figure 15 shows the 4 form of the transfer function for the two -hole mask. We have used 
square holes in this example to simplify the evaluation of H^(k). 
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Pig. 13. Lens geometry. 

The field over the focal plane of an aberration-free lens is closely approximated by 
the relationship 

i f* -j2Tru • x/XL 

. E v (u,t) = xl J dxw(x) E v (x,t) e . , (84) 

where w(x) is the transmittance function for the lens and may be possibly complex (phase- 
shading). The intensity of the focal plane field is 

Av f C * -j2iru- (x-x')/XL 

I (u. t) = y \ dx \ dx' w(x) w(x') T (x,x',t) e , (85) 

v (XL)^ J J v 

where T v is the coherence function measured per unit optical bandwidth, and the 
factor Av accounts for the integration across the correlation bandwidth. 

The spatial Fourier* transform of this image (with a sign reversal employed to 
account for the image inversion), is simply 

+j2irk-u/L 

«^ v _(k. t) = du I(u,t) e 

= Av y dx w(x) w(x-Xk) r y (x, x-Xk,t). (86) 

In the normal use of a lens the integration on x extends over the full aperture. We 
can localize the measurements of coherence by an appropriate choice of transmittance 
function w(x). Suppose, for example, that we mask the lens except for two small holes, 


33 



Besides localizing the coherence process in the x variable to a single correlation 
area, a receiver should also integrate r over a range of displacement values. This 
operation is mathematic ally ~ equivalent to restricting the angular field of view to just 
the region of the object, to exclude as much of the background as possible. We do this 
in the focal plane by inserting a field stop, another mask with a small hole matched to 
the image. 

When a telescope is used in this manner it is commonly referred to as a Michelson 
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stellar interferometer {see Born and Wolf ) which is used in the measurement of star 
diameters. The spacing between the openings in the mask is increased until the fringes 
disappear. At this point the wave -front coherence has dropped, thereby indicating that 
we have passed the uppermost significant spatial frequencies in the object, and from 
such knowledge we can deduce the approximate dimensions of the object. 

Spatial frequencies greater than those admitted by the aperture of the lens can be 
sampled by using mirrors to extend the baseline (see Fig. 16). 

! . j 

\j- — ^ ^ ^ 



Pig. 16. Baseline extension. 

Changing the spacing of the outermost mirrors is a convenient way to sweep across 
the spatial-frequency spectrum, since the fringe spacing remains constant in the focal 
plane. 

Finally, we might say a word regarding the practical significance of the underspread 
condition (section 2. 3. 3). When the holes in the mask have the characteristic dimension 
Ax, the resultant diffraction effects spread the image in the focal plane to an apparent 
angular width of approximately X/Ax. That is, neither the amplitude nor phase of the 
fringe pattern can change on a scale finer than this. Thus if we are to resolve the angu- 
lar correlation region, we require X/Ax < Az. 

2. 5 NOISE IN THE MEASUREMENT OF MUTUAL COHERENCE 

We have just considered various interferometer schemes for sampling mutual coher- 
ence. We want to sample the mutual -coherence process in all of the relevant variables — 
time, optical frequency, aperture space and displacement — in order to retain all of the 
statistically relevant information contained in the incident field. Such a sampling should 
give us maximum use of the spatial and temporal channel diversity. We have considered 
several schemes for sampling the coherence function. Our object now is to study 
the additional disturbances that the quantum nature of radiation introduces into the 
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interferometer measurements. We shall have then completed the channel modelling, 
and can go on to the performance of optimum processing schemes for parameter 
estimation. 

Let us look at the statistical behavior of the measurement noise that arises because 
the basic principle of energy detection is a counting of photons whose arrival is stochas- 
tic in nature. If the intensity of the radiation is very low and we employ a detec- 
tor with a large bandwidth, then individual photons (or more properly photoelectrons) 
can be resolved. Usually, however, the detector current represents an integration 
of a large number of such current pulses. Let us suppose that the detector has 
a minimum integration or resolving time, r. We now ask for the statistical behavior 
of the detector output, still conditioned on knowledge of the intensity of the object 
and the amplitude and phase disturbances of the atmosphere (a,b). The detectors 
should be chosen to have a resolving time less than the correlation time of the tur- 
bulence. 

There are two random phenomena that yield uncertainty in the detector output. The 
first is the random photon- detection process, the second is the Gaussian nature of the 
incident radiation. 

Let us subdivide our resolution period r into small intervals 6^, 6^ ... 6^ (N=t/ 5), 
where 6 is shorter than the coherence time of the incident field. If we condition on 
knowledge of this field, the corresponding counts c^, . . . c^ are independent, Poisson- 

distributed, random variables with a mean and variance 

c. = (c.-c.) 2 = g 6 Au 2 |e.| 2 , (89) 

where 


2 

Au is the area of the detector 

|E.| 2 is the instantaneous power flux density of the incident field 


and 


q (quantum efficiency) 
g = — 

vh (frequency X Planck's constant), 


where g is in counts/joule. 


The total count in the period t then is 


N 



i=l 

% 


Averaging over the Gaussian statistics gives a mean count of 


(90) 


( 91 ) 
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( 92 ) 


c = g6Au 


N 

l |E i |: 


1=1 


A ^ 

= gTAu I, 

where I is the average power flux density (the double bar is used to indicate averages 
over the two random phenomena). This average count then represents the signal com- 
ponent in the output. The variation from this average count we call the noise, and it 
has variance 

= N N 

(c-c) 2 = g 2 6 2 Au 4 ^ ^ |E(t.)i:<t.) 1 2 + g T Au 2 I. (93) 

i=l j=l 

Under the assumption that the resolving period t is many times longer than the corre- 
lation time of the incident process, one of the summations in (93) can be extended to -H*> 
without significant error, and we obtain 


, =,2 2 A 4 

(c-c) - g tAu 


p+00 

^0 


r(t) | 2 at + gTAu 2 i, 


(94) 


where F(t) is the stationary temporal coherence function for the field. We can define a 
coherence time r for the process as 


Jj 00 [r(t)j 2 dt 


where 


lno)| 2 . 


(95) 


We thus obtain 


. =.2 . 2 t , 2 2. 4 t 2. / . 

(c-c) = grAu I + g t Au I (t /t), 

c 


(96) 


where 


c 


gT Au 2 I. 


(97) 


We see that the noise variance comes from two independent effects. There is first 

2 

the shot-noise term gTAu I. It would have the same value if the field were constant, 
rather than random. The second term, g t Au I (t /t), results from the residual 
uncertainty in the sample average intensity when approximately t/ t independent sam- 
pies are considered. Generally, t/ is extremely large (easily of the order of a mil- 
lion). Thus this 11 phot on -bunching” term is usually insignificant in comparison with the 
shot-noise term. The effect is really the same phenomenon as that in the classical 
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approach when we considered the variance between the sample-average and ensemble- 
average coherence functions. 

Let us now examine the correlation of the noise for a pair of observations. Consider 
a pair of fields Ej(t) and E^(t) related by a temporal cross correlation function T ^ 

Each of these fields is detected on separate counters giving counts 


c = Cl + c 2 + 


"N 


C 1 = c'j + + 


pi 

N 1 


( 98 ) 


where the resolving period has again been subdivided. We want to consider the correla- 
tion 


(c-cHc'-c 1 ) = 


l KJ - 

i=l 


gT Au^Ij 


N 


l [ 1 ] - 

L3=l * 


gTAu 2 I 2 


N N 


i=l j=l 


( 99 ) 


Now, again extending one pair of the summation limits to +<*> and defining a correla- 
tion time between the two processes as 


T - ' 

C 




( 100 ) 


we obtain 


(c-c)(c'-c') = g^uVjr^ 2 ( 0 ) f 2 (t c /t). ( 101 ) 

This crosscorrelation term for the noises is purely a 11 photon-bunching 11 effect, shot- 
noise components being independent. When there is no coherence between the fields the 
correlation between the noises is strictly zero. This is the case when the two time 
intervals are disjoint and each is much larger than the coherence time of the radiation. 

It matters not whether we consider two different detectors or the same detector. 

When we have field coherence there is some noise correlation but, since this corre- 
lation is of the same order as the “photon-bunching 11 term in the variance function, it 
is generally quite insignificant. 

We can, therefore, for the most part of the succeeding analysis, disregard all but 
the shot-noise term. The other term does have significance in the Hanbury Brown and 
Twiss intensity interferometer, which will be briefly considered in Section III. With this 
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technique, field coherence can be deduced from the correlation of intensities, because 
of the "photon-bunching 11 term* This, however, would require a very large detector 
bandwidth in order to be practical. 

The fact that the dominant shot-noise components are uncorrelated for disjoint obser- 
vation times and disjoint detector areas, simplifies considerably the analysis of the noise 
in our coherence measurements. 

Consider the result of correlating an intensity fringe pattern I(u) incident on an array 
of detectors with spatial* frequency cos (2irk*u). The result of this operation we denote 






I(u) cos (2irk • u) du. 


(102) 


We determine S as 
c 

S c = ^ cos (2Trk*m), (103) 

i 

where we have denoted by c^ the photoelectron count from the detector element at 



point m (see Fig. 17). 

From Fq. 92, it follows that 


= 1° i 


cos (2'irk’m) 


= gAT ^ Kuj) COS (2Trk‘U^) Au^, 


(104) 
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where the integration time is now AT, the correlation time of the atmospheric distur- 
bances. We represent the averaging by a single overbar, as it is no longer necessary 
to distinguish the two imbedded random processes (the Gaussian radiation and the 
Poisson photodetection process). The summation in (104) can be replaced by an inte- 
gration (under the assumption that the detector elements are small), and we obtain 

§ c = gAT ^ I(u) cos (2irk * u) du. (105) 

In terms of J?( k), the spatial Fourier transform of the pattern, becomes 

S = gAT Re [/(k)]. (106) . 


Likewise, when we correlate against sin (27rk*u) to obtain the quadrature component 


S cc 
s 


I 


I(u) sin 2irk • u) du. 


(107) 


we obtain 

S g = gAT Im [/(k)]. (108) 

The effect of the shot noise resulting from the random emission of photo electrons 
gives two zero-mean corresponding noise variables n^ and n g . The mean- square value 
of each of these is the same, and is obtained as 


2 ^2 
n^ = (s -S ) 
c c c 


■n 


(c.-c.)(c -c.) cos (2irk/u.) cos (27rk-u.). 
i 3 3 ^ 3 


(109) 


From Eq. 101, we found that the covariance of the counts on disjoint detectors is due 
entirely to the "photon bunching," that is, to the statistical structure of the received 
radiation, and is negligible for the relatively long integration times employed. Thus 
(109) reduces to 


2 \ , - .2 2 /n , 
n = / (c.-c.) cos (2Trk-u.) 

c Av 1 1 “ “i 

i 

which, with the aid of (96) (retaining only the shot term), becomes 
n 2 = gAT ^ I(u) cos 2 (2Trk*u) du. 


(HO) 


( 111 ) 
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Similarly, the noise on the quadrature component is 


= gAT y I(u) sin^ (2irk*u) du. 


(112) 


These can be rewritten 


n 


gAT 

2 


[/(Q) + Re jP (2k)] 


and 


n 


2 

s 


gAT 

~ir~ 


[ JP (O) - Re </(2k)]. 


(113) 


(114) 


In the case of the simple interferometer forms (sections 2. 4. 1 and 2. 4. 2) the fringe 
pattern does not contain any component at the double frequency 2k. Moreover, in the 
case of a lens with a two-hole mask, the double frequency is usually beyond the limit 
of the MTF (Fig. 15). Thus, we can drop the second term in (113) and (114), which 
leaves 


2 

n 

c 



gAT 

2 


JP( O). 


(115) 


When the double- frequency component is absent, it also turns out that n and n are 

c s 

uncorrelated. 

We now apply the results of Eqs. 106, 108, and 115 to the simple interferometers or 

the masked lens. The squared magnitude of the mean complex signal S = S + jS is 

c s 


|s | 2 = g 2 AT 2 Av 2 Ax 4 l® v (k) I 


2[a v (x, t)+a v (x-Vk, t)] 


(116) 


Similarly, the mean- square value of the magnitude of the complex noise compo- 
nent n is 


n 


2 = gAT Av Ax 2 


© 

v 


(O) e 


2a v (x, t) 


N v « 


'][« 


2a (x-\kt) 

© e v “ + 

V 


N v « 




(117) 


The variance of the shot-noise variable is dependent upon both the background -radia- 
tion N^, the object radiation and the atmospheric fading a v (x, t). When we come to 
study performance limitations, we shall simplify the analysis by supposing that is 
sufficiently large that it is the dominant component in the noise variance. This then 
removes the dependence of the noise variance on the atmospheric process and- the object. 
Furthermore, we assume in Section III that the noise variable n represents an 
average over a sufficient number of counts that a Gaussian approximation is valid. 
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The signal-to-noise ratio then is 
| S | 2 gAT Av Ax 2 

|© v (k) i 2 ex P { 2 t a v (x> t)]}. (118) 

This signal-to-noise ratio was calculated for a scheme in which a spatial fringe pat- 
tern is formed. If we use a moving mirror to obtain a periodic detector signal, the 
signal-to-noise ratio in the coherence measurement will be exactly the same. Scanning 
the fringes with the ruling device of section 2. 4. 1 does result in a slight loss in signal- 
to-noise ratio. Some of the light energy is blocked by the opaque regions, which results 
in a loss of a little more than a factor of two in signal-to-noise ratio. 

2. 6 COMMENTS 

The results of section 2. 5 give the signal-to-noise ratio for the measurement of the 
mutual-coherence process, isolated to a single correlation time AT, correlation band- 

z 

.width Av, pair of correlation areas Ax and field of view ft. In order to retain all 
information in the received phenomenon (at least in the classical sense), we should 
simultaneously sample the mutual- coherence process for all such available frequency 
bands, time intervals, pairs of correlation areas, and spatial frequencies* Conceptually 
at least, simultaneous sampling in frequency could be done by first' diverting the light 
into a number of channels, one for each frequency band, with the use of some system of 
lossless interference filters. There would be no loss in the signal-to-noise ratio given 
by Eq. 118. There is no trouble in sampling all time intervals. 

.It is in sampling over aperture position and displacement (spatial frequency} that 
trouble arises. First, however, consider the displacement as fixed. The field incident 
upon a particular correlation area should then be interfered only with the fields through 
two other such areas. If this is done by using different interferometers, the signal-to- 
noise ratio is then reduced by a factor of two, because of the necessary energy division. 

This is not too serious; however, when we try to sample over the full range of spa- 
tial frequencies we encounter problems. Now we must interfere the field on a particular 
correlation area with those on all other areas. If we have M areas, we might use beam 
splitters to funnel the energy to M-l additional interferometers. Unfortunately, the 
light energy is then split M ways, with a subsequent M-fold decrease in signal-to-noise 
ratio on each interferometer measurement. 

The problem is more subtle than this, in that we can devise schemes of introducing 
a spatially varying frequency shift to the incident field, using systems of moving mir- 
rors, etc. There is no division of the light energy, and indeed only a single detector 
need be> employed. The coherence at various values of both position and displacement 
show up as orthogonal frequencies in the spectrum of the detector output. In terms of 
our noise model (section 2 . 3 ), this means that the signal component of the detector 
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output is unchanged from a system using only a single pair of correlation areas. Unfor- 
tunately, the noise takes an M-fold increase because the total energy incident on the 
detector is now M times as large. Thus, no matter how ingenious we try to be, it 
appears that we are faced with the fundamental problem of compromising between 
a small set of spatial frequency samples with minimum noise and a large set with much 
higher noise. 

In Section III we shall consider the limitations on the best possible strategies of esti- 
mating the phase and amplitude of the Fourier transform of the object pattern, using the 
data from a particular set of measurements. After understanding what can be done with 
a particular set of measurements, we can determine the most reasonable compromises 
for selecting a sampling scheme. 
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Ill* PERFORMANCE LIMITATIONS AND SIGNAL- 
PROCESSING TECHNIQUES 

We have determined the statistically sufficient set of measurements that a receiver 
should perform upon the atmospherically distorted incident field to retain all relevant 
information for optimal reconstruction of the image. We shall now explore the manner 
in which these measured values should be processed to form the image estimate. First, 
we determine the best possible performance by means of the Cramer -Rao (C-R) inequal- 
ity which lower-bounds the variance of any unbiased estimate. Second, we demonstrate 
receiver structures that either attain this lower limit or approach it asymptotically. 
Finally, we discuss various practical realizations of these structures. 

We have shown that sufficient field measurements are obtained by a bank of parallel 
interferometers that measure the mutual coherence process, which is defined as 

r v (x, x+y, t) = <E v (x,t)E v (x+y,t)}. (119) 

The sample averaging denoted by angular brackets is a local integration, or smoothing 
of the mutual coherence function, where 

v the optical frequency is integrated over the optical correlation 
bandwidth Av 

t the time is integrated over the correlation period At 

2 

1 x the aperture coordinate is integrated over the correlation region (Ax) 

2 

y the displacement variable is integrated over a correlation region (Ay) , 
equivalently 1 a restriction of the field of view to a solid angular 
region Q = (X/Ay) . 

These correlation regions, for the most part, are the corresponding correlation 
regions for the phase and amplitude disturbances induced by the turbulent atmosphere. 
Because of the smoothing, the mutual -coherence process can be sampled in v, t, x, and 
y (with a sampling grain equal to the forementioned correlation regions), rather than 
being retained as a continuous function. 

Under the isoplanatic condition, these samples are related to the object, the atmo- 
spheric disturbances, and the shot noise as 

? v (x, x+y, t) = ATAv{Ax) 2 © V (y/M exp{a v (x, t)+a v (x+y, t)+j[b v (x, t)-b v (x+y, t)]} 

+ N v ATAv(Ax) 2 S2 60,y + n y t , (119a) 

where 

© v (k) is the spatial Fourier transform of the intensity pattern of the object 
about optical frequency v 
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<v. 

N is the background radiation level 

a^(x, t) is the atmospherically induced log amplitude disturbance 

t> v (x, t) is the atmospherically induced phase disturbance 

n . is the complex shot-noise variable, 
v, x, x, y 

We assume that a and b are jointly Gaussian processes, and that the n . are 

v , y 

independent complex Gaussian-noise variables. For ease of analysis, we suppose that 
the total received energy {mutual coherence process at zero displacement) is dominated 
by the background radiation. The variance of the shot-noise variables, which depends 
upon the total received energy, is then essentially independent of both the object and the 
atmospheric disturbances. 

An estimate of the object may be either a full restoration of the intensity pattern or 
simply an estimate of a small set of parameters, perhaps the outline shape or dimen- 
sions of the object. In the case of a full restoration we need to know all of the complex 
values of © v (k), whereas in the case of parameter estimates a few values may suffice. 
Indeed, we are often satisfied. with only the magnitude of the transform at a few points. 
For example, in estimating the diameter of an object, one is often satisfied to determine 
the spatial frequency at which the magnitude of the spatial transform drops to some 
prescribed fraction of the .zero -frequency component. 

It is useful to separate the possible estimation problems into two categories, those 
that require only the amplitude of the spatial Fourier transform, and those that require 
both the amplitude and the phase. 

Consider the estimation of the magnitude or phase of a particular © v (k), and let us 
look at just the samples of T for displacement y = \k. Furthermore, assume that these 
samples extend only over time. We shall now find the lower bounds to the variance for 
any estimate, and then demonstrate the receiver structures that approach these limits. 
With this understanding of the elementary form of the problem, we then consider the gen- 
eral case for which we may have available extra information in additional aperture coor- 
dinate, displacement, and optical-frequency samples. 

We make one further approximation in order to shorten the problem, and treat the 
successive samples as being statistically independent. For the shot-noise components 
this is certainly valid, as the observation times are disjoint. It would be presumptuous, 
however, to suppose that atmospheric disturbances are independent, just because the 
samples are separated by the correlation time AT. In general we expect some smoothly 
decreasing degree of dependence as the separation between samples is increased. 

Nevertheless, from our experience we believe that the results obtained under the 
independence assumption will usually reflect the behavior of the channel with sufficient 
accuracy for engineering considerations. The effect of the actual intersample depen- 
dence can usually be taken into account with 
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1. An effective signal -to -noise ratio which may be somewhat greater than that of 
each sample. 

2. An effective number of independent samples that may be somewhat less than the 
value that we shall use, N = T/AT (T being the total duration of the observation). 

For convenience, we denote the successive time samples of T as 

r i = ryx, x+Xk, tj) 


= a ! / 2 


e^e 




n i? 


where 

A 1/2 =|® v «|e v 

cp = phase angle of © v (k) 

~ a v <5» + a v <x+Xk, t i ) - 2a v 

^ = b v (x, t ± ) - b v (x+\k, t i ) 

n i - n v, tj, x, x+Xk' 


( 120 ) 


( 121 ) 


Let us normalize the noise samples so that each component (real and imaginary) has 
unit variance 


E[Re<n.)] 2 = E[lm (m)] 2 = 1. 


( 122 ) 


The parameter A is thus proportional to the squared magnitude of the spatial Fourier 
transform component. With the normalization of the noise, A is a measure of the aver- 
age signal-to -noise ratio in the output samples. 

In the succeeding analysis we frequently encounter the ensemble average values of 
na, jnb. 

e and e , which we denote as D and B (n an integer). Under the Gaussian 

n n 

assumption 


n 2< r 2 /2 

»n=* na =e a 


B 


n 




(123) 


2 2 

where o* a and are the mean-square values of the log amplitude and phase disturbances. 
The mean values of a v (x, t) were subtracted in the definition of a (Eq, 121), and absorbed 
into the definition of A in order that have zero mean. 
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3.1 LOWER BOUND TO THE VARIANCE OF ANY UNBIASED 
AMPLITUDE ESTIMATE 


The Cramer.-Rao inequality gives as the smallest possible variance of any unbiased 
~ I 


a 13 

estimate A. 


Variance A ^ 


E 


9 


8A ln P( r !» r 2 ' • - r ; 


I 2 ’ 

4 1 ) 


(124) 


or equivalently 


Variance A ^ 


-1 


r 2 n * 

9 ^ 

. — 2 ln p(r l’ r 2* ■ • r N/ A ' ^ 
_3A 


(125) 


E 


where the expectation (denoted by the operator E r ) extends over the random variables r^. 
The spatial phase <f> is regarded as a known parameter. With N independent samples, 
this inequality becomes 

^ 1 


Variance A > 


N 
S E 
i=l 


r[A ln P(r i /A ’ «] 


or 


-I 


Variance A 5* 


N p 2 ”1 * 

s B B-ln p(r./A, <{>) 

i=i r LaA z 1 J 


( 126 ) 


(127) 


The problem encountered in directly evaluating the C r R inequality is the complicated 
form of the conditional probability resulting from the mixture of normal and log-normal 
forms. We are forced to take a piecemeal approach to the evaluation, using different 
approximating techniques over various ranges of the parameter A. 

Consider, first, the two asymptotic cases of very large, and very small values of A. 

3. 1. 1 Very Strong Signal 

When A is very large, we completely disregard the shot-noise component, thereby 
making each r. a complex log-normal random variable. Taking the logarithm of each 
sample then reduces the problem to the familiar one of constant signal in added Gaussian 
noise. Application of the C-R inequality then yields 

„ • £ 4 <r 2 

Variance A ^ a / 140 i 

— jz — *ir • (128) 


3.1.2 Very Weak Signal 

The opposite extreme of very poor signal-to-noise ratio is harder to work out. In 
Appendix B we derive an asymptotic expansion of the right-hand side of the inequality 
of (126). Under the condition that the phase variance is large (B « 1), we obtain the 
small- signal approximation 
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where only the first two significant terms of the power-series expansion of the denom- 
inator have been retained. 

If we suppose that N (the number of independent observations) is the same in all 
cases, we can write 

Variance A ^ 

A 2 " N ’ 

where the function G{A) provides a convenient means of comparing the performance lim- 
itations. The general form of G(A) for the two extreme ranges of A is shown in Fig. 18 
(log-log scales are employed). 

The first result (inequality (128)) is valid for all values of A, but is not useful except 
for A large, as it reflects only the limitation caused by the turbulence, and ignores the 
effect of shot noise. 

The other extreme case (inequality (129)), shows a variance which decreases 
inversely as A in the limit of small A. For larger values of A such that A >D^B^/D^ 
the function changes to an inverse quadratic in A. It turns out that this is the essential 
behavior over a very large range of A, but does not necessarily follow from (128) 
because we simply ignored the higher order terms in the expansion. 

3. 1. 3 Intermediate Signal 

We fill in the intermediate region of the bound in two more steps. 
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a. Uniform Phase and Constant Amplitude Approximation 


The continuation of the low -signal form (129) to higher values of A is obtained by 
the following approximation. First, we suppose that there is no fading and replace the 
a i 1 /2 

term e in Eq. 120 by the factor m This ensures that the same mean-square value 

of the signal component is retained. Furthermore, we suppose that the phase probability 
density is uniform over the closed interval 0 to 2tt, which would seem to be reasonable 
when the variance of the phase is large. The C-R inequality can then be approximated 
as (see Appendix B). 


Variance 


A 

A 


16 


A 


N SD^+SD^+D^A 4 


(130) 


Only when the signal-to-noise ratio is quite large (A> 1 /D^) do the third- and fourth- 
order denominator terms give a significant contribution. Thus over a large range of A, 
the inequality (130) can be reduced simply to 


Variance A 


nd|a 2 


(131) 


This inverse quadratic behavior is almost the same as that given by the second term of 
(129) (except for a factor of two). Indeed, we expect that (129) is the more accurate, 
since one of the approximations for (130) was the elimination of the fading. 

The implications of the uniform-phase approximation must also be considered. We 
note that the expression (130) does not tend to a l/A form when the signal is very weak, 
as predicted by the expansion of (129). The uniform-phase approximation is the root of 
this anomaly. It comes about because the mean component of each r^ is destroyed under 
a uniform-phase approximation. As we shall see eventually, the optimum receiver in 
the limit of very small values of A is one that treats this small mean value as a known 
signal in an added, zero -mean background, and ignores any information contained in the 
rest of the sample. 

Admittedly, the uniform -phase approximation does not give a rigorous lower bound 
to the problem. It is our opinion, however, that the true Gaussian nature of the phase 
has an appreciable effect only when the signal is very weak and when the very small 
mean component in the signal can be exploited. 

For larger values of A the cubic and fourth-order terms become significant, and 
the utility of the bound is reduced. Rather than worry too much about these terms, we 
turn instead to the last step in the bounding procedure. 

b. Known Phase and Amplitude Approximation 

A strict lower bound can also be obtained if we assume that both a.. and b^ are 
known to the receiver a priori but are still random. We then have a case of signal 
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in added Gaussian noise. The C-R inequality is easily evaluated to give 


Variance A 
A 2 


N 

A 2 e 
i=l 


2a. 

i 


which when' averaged over yields 

Variance A 4 
a 2 - HAD., ' 


( 132 ) 


where we have exploited the well-known result that for a positive random variable x, 
(l/x) » (l/x). 

The four pieces that approximate the C-R lower bound are summarized in Table 1. 
The form of this bound is sketched in Fig. 19. We divide the A axis into four regions. 
Beginning in region 1 where the signal level is very high, the performance is indepen- 
dent of the signal- to- noise ratio A. When A is somewhat smaller, the variance varies 
inversely as A (region 2). At still lower values (region 3) the behavior is an inverse 
quadratic. Finally, in region 4 the signal is extremely small, and the variance follows 
inversely as A. 

3.2 RECEIVER STRUCTURES FOR AMPLITUDE ESTIMATION 


Let us now look at various receiver structures that we can use in order to' actually 
achieve performance either as good as or almost as good as the lower limits specified 
by the C-R inequality. We shall assume that the receiver has enough a priori knowledge 
regarding A to know the general region of operation (strong signal, intermediate signal 
or low signal). 

3. 2. 1 Fringe-Tracking Receiver for Strong Signals 


The first receiver to consider is the one that tracks the fringes in the interferometer. 
That is, it is able to obtain accurate measurements of the spatial amplitude and phase 
of the fringes in each coherence interval AT, without shot noise significantly distorting 
the result. 

The log amplitude and phase of the complex samples are conveniently expressed 
as the two complex components of In r^. 

The maximum- likelihood estimate for A follows immediately as 



(133) 


This maximum-likelihood estimate for A has a bias 



- A = A 



(134) 
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When N is large, the bias becomes negligible. Moreover, we can remove the bias com- 

( 2 ^ ) 

pletely simply by dividing the estimate by \^e - 1/ , but the estimate is no longer 

maximum -likelihood . 

The fractional mean- square error in the maximum -likelihood estimate is 


[Vr A l 2 


A 



+ e 


8c 2 /n' 
a' 


(135) 


When N is large enough that the error is small, 10% or less, the exponentials can 
be expanded to first order, and we have 




(136) 


This is identical to the C-K lower limit (128), and thus this amplitude estimate is said 
to be asymptotically efficient with N. 

With noise present, the performance of the fringe- tracking receiver is naturally 
degraded. As long as the noise is small, we can approximate 


In 


1 +• 


n 


L A l/2 e a+j(b+<t>)J 


nA'l/ 2 e - a -J‘ (b+4,) . 


(137) 


The phase term e - ^ b+< ^ only rotates the noise coordinates, and can be ignored, as far 
as the performance goes. 

Each complex component of In r^ then has an additional random-noise variable 
A 1 / 2 e &1 Re (n.) 


and 


■ 1/2 


-a. 


A A/ ^ e 1 Im (n.). 


(138) 


These are uncorrelated with the phase and amplitude disturbances, and provided the 
error in the estimate is small (N sufficiently large), the resulting error is 


[jU] 2 *['» +D 2/ A 


N 


(139) 


The linear approximation to the logarithm holds only when n. is quite small, ‘Indeed, if 
y a. 1 

|n^| 2* A 1 ' e 1 , the second-order term exceeds the first-order term in the series 
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expansion (137), and the linear approximation is not justified. 

The nonlinear noise effects are not nearly as disastrous as in the phase estimation 
problems, but a bias will be introduced into the estimate, and the performance will not 
be as good as that for the next receiver structure. 

3. 2. 2 Square-Law Receiver 

If we sample -average the squared magnitudes of the samples r^, the ensemble aver- 
age of the result is 

N 

H = ^ ^ |rj 2 = AD, + 2, (140) 

i=l 


where the factor 2 is a bias arising from the noise energy. An estimate 


A 

A = 


1 1 


N 


- 2 


(141) 


is then unbiased, and has a fractional mean-square error of 


[A-A] 2 


N 


«-» 1 T-'l A * 


D 2 A ‘ D 2 A 2 j 


(142) 


When the signal is very strong, A > 4(p^-D 2 )/D 2 . The error is independent of the 
signal strength and is approximately 



(143) 


2 

The corresponding C-R limit is 4cr /N. Thus the square-law receiver is in gen- 

a 

eral not efficient whenever the signal level is high (region 1). In the case of very 
slight fading small) a first-order expansion of the exponential suffices, and we 
find that 


[A-A] 2 4cf 2 

3> 

A 2 " N ' 


(144) 


which is the same as that given by the C-R inequality in region 1. 

When the signal is somewhat weaker, the inverse linear term dominates the error 
expression, and we have 
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(145) 


[A- A] ‘ 
A 2 


“ « [“ J 


D 2 < A < 4; 


D 4 - D 2 


D 


2 J 


This is identical to the limit given by the inequality (132), and the square -law 
receiver is thus efficient in region 2, 

When A is very small, the fractional mean-square error behaves as 


[A-A]' 


N 




A < D' 


-1 


(146) 


This can be compared with approximations (130) and (132) for the C-R limits in region 3 
4 


and 


N 


1_ 

N 


A Z Dl 


(second term of asymptotic expansion) 


a 2 d. 


(a constant, b uniform, 0 to 2-rr). 


If we suppose that the first limit (130) is the more accurate, we find that this square- 
law receiver can actually do as well as could ever by expected in region 3, 


3,2.3 Cross correlation Receiver 


We can make a slight modification to the square-law receiver which avoids the prob- 
lem of noise bias. The bias must be subtracted from the sample average of the squared 

/% 

magnitude of the samples before we obtain the estimate A. This number would either 
be known a priori or would need to be estimated, but by going to a cross correlation form 
of the receiver we can eliminate the average noise energy in the output. Suppose that 
the fringes formed in the interferometer are measured by two detectors. We then get 
pairs of samples which have an identical signal component but independent shot-noise 
variables. In a crosscorrelation of these samples, the noise products then tend to zero. 
Because of the energy division, the signal levels and the variance of the shot noise are 
half the value for a single detector. Normalizing to unit variance of the noise compo- 
nents yields sample pairs of the form 


r i = 


1/2 a i + i (b i +< » 

i 1 e 

42 


+ n. 



(147) 
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The ensemble average of the sample crosscorrelation product is 


N 



i=l 


AD, 


and thus an unbiased estimate can be obtained as 


A 

A 


2 1 
d 2 n 


l 

i=l 


❖ 

r i v r 


( 148 ) 


(149) 


The fractional mean-square error in such an estimate. is 




8 

ad 2 


16 




( 150 ) 


As should be expected, when the signal level is very high the variance in the esti- 
mate is the same as with the square -law receiver. In the two other regions the behavior 
is functionally the same as for the square -law receiver, but we suffer a penalty, in that 
A must be twice as large in order to achieve the same results. That is, we forfeit a 
factor of 2 (3 dB) in the effective incident light flux. Nevertheless, when the signal -to - 
noise ratio is poor one is vulnerable to a gross error with the square -law receiver, 
unless the noise bias is known very closely. 


3. 2. 4 Signal Integrator 

In order to attain the very small signal performance predicted by the C-R inequality 
of (129), we must employ one final receiver form. This is a receiver that considers the 
residual mean component in each sample r^ as a known signal imbedded in a zero -mean 
random disturbance. If we form a sample average of the r^, we obtain a random vari- 
able R that has an ensemble -average value 


N 



i=l 


= A 1 / 2 DjBj. 


-A 

An estimate A can then be formed as 


(151) 


A 

A = 


D i B i 


*1 

7 


2 


r i 


The ensemble average of this estimate is 


( 152 ) 
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A 

A = A 


2 


(153) 
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2 
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N 

L b x J 



ND 2 B 2 


A bias exists in the estimate, but it tends to zero with increasing N. Moreover, it 

2 , 2 . 

can be removed by first subtracting the noise contribution z/ND^B^, and then dividing 
by 1 + (1 /N)(d.,/B 2 D 2 -i). 

The evaluation of the mean-square error, though straightforward, is extremely 
tedious. After going through considerable algebraic manipulations, and elimination of 
higher order terms of l/N which are negligible when the final error is small, we obtain 


[A-A] 2 

A 2 


N 


AD 2 B 2 


2P 2 

D f B l 


+ 2 Z 4 


(154) 


Moreover, in cases of practical interest, B^ and B,, are both very small, and the 
expression reduces to 


[1-A] 2 

A 2 


_1_ 

N 


AD i B ? 


2D Z 1 


(155) 


For very small values of A the second term is negligible, and 


[A-A]' 


1 

N 


A»i B i 



(156) 


which is identical to the asymptotic result for the C-R inequality (region 4). 

2 Z Z 

Throughout the range A > D^B^/D^ the integrating receiver is suboptimum, and 
indeed when A > 2/D,, the performance tends to saturate at 




2D, 


D i B ? 


(157) 


This is far greater than the saturation values for both the fringe -tracking receiver 
{Eq. 136) and the square-law or crosscorrelator receiver <Eq. 143)* The integrator 
receiver considers only the small mean component as useful signal, and regards all the 
rest as noise. Even when the shot noise is absent, the tremendous difference between 
this mean component and the rest of the signal creates a very small effective 11 signal - 
to -noise ratio. 11 

Although the integrating receiver is inferior at the high values of A, it is interesting 
to note that as the signal becomes weaker, we eventually reach a point where the linear 
behavior of the variance with respect to A triumphs over the noise -enhancement 
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A 

Table 1. Cramer-Rao lower bound to variance of A. 


A 

Table 2. Receiver characteristics for A. 


RECEIVER CHARACTERISTICS 
[£ - a] 2 /A 2 « (1/N)G(A) 

Gl ( A) = 4 [ + D 2 /A] 

FRINGE 

TRACKER 

weak noise 

p 

~ D[i-D2 +4 + 4 

G (A) = 

d| ad 2 a%| 

SQUARE 

LAW 

noise power 
known aprlori 

Dk-dI + 8 +16 

gL (a) = — 

d| ad 2 A^l 

CROSS- 

CORRELATOR 

no knowledge 
of noise 

- , ^ 2 D 5 

. Gt, (A) - + 

AD^bJ dJ Bf 

ii li 

INTEGRATOR 


CRAMER - RAO 101® BOUND 
Variance A/A^ ^ (l/N) G(A) 

REGION 1 
noise = 0 
strict 
lower bound 
for all A 

g 1 (a) = Mr a 2 

REGION 2 
a,b, known 
strict 
lower bound 
for all A 

g 2 (a ) = 

“2 

REGION 3 
a constant 
b uniform 
approximation 
to bound 

16 

G,<A) ~ 

8d |a 2 +8d |a 3 +d^a^ 

REGION k 
signal very- 
weak 

asymptotic 

approximation 

h 

V A) =-5-£ z-2 

Df B| A + D2 A^ 





characteristic of the square -law operation. 

The performance results for the three receiver forms are summarized in Table 2, 
and can be compared with the Cramer-Rao lower limits of Table 1. The general func- 
tional forms of the various error expressions are plotted in Fig. 20. 

3.3 LOWER BOUND TO THE VARIANCE OF ANY UNBIASED 
PHASE ESTIMATE 

In an almost identical manner, we can derive a lower bound to the variance in the 
phase estimate <j>. 

3.3.1 Strong Signal 

When the signal is very strong, 
lower bound in region 1. 

2 

a a b 

Variance <f> > 

3.3.2 Very Weak Signal 

At the other extreme of very weak signal (region 4) we again make an expansion of 
the right-hand side of the C-R inequality where the parameter A is now considered to 
be known. Thus the asymptotic form of the lower bound becomes (Appendix B, sec- 
tion B. 3) 

Variance <i> 5* (159) 

NADj Bj 

3.3.3 Intermediate Signal Level: Known Channel Phase and 
Amplitude Approximation 

The intermediate region 2 can likewise be evaluated in a manner similar to that used 
in the amplitude -estimation problem. We assume that a and b are known to the 
receiver. Conditional on this knowledge, the C-R inequality is simply 

A 1 

Variance <j> > (160) 

S A e 1 
i=l 

An average over the a^ then yields 
A 1 

Variance <f> > NAD . (161) 

The intermediate behavior in region 3 is the most difficult to obtain. We strongly 


we can discard the noise terms and thus obtain the 

(158) 
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suspect that when the phase variance is large there is a severe threshold effect. That 
is, there seems to be little intermediate region. Provided the signal dominates the 
noise, a relatively good phase estimate is possible.* As soon as the noise is dominant, 
however, only the residual mean component of the signal is of use in phase estimation, 
and we operate in region 4 with the much larger variance given by (159). 

3. 3.4 Noisy Phase Reference and Constant Amplitude 


We can derive a lower bound to the variance in the estimate of the phase difference 
between two signals that suffer the same atmospheric disturbance but have independent 
added noises. This bound will be useful later when we consider a procedure that forms 
a phase estimate by sequentially estimating the phase difference between pairs of neigh- 
boring spatial frequencies. 

Let us suppose that in addition to the samples r. we are also given a set of noisy 
phase reference samples v. of the form 


v i = A 


1/2 


+ im, 


( 162 ) 


where m. is an independent noise voltage having the same mean-square value as n.. The 
phase <f> 'is then the phase that would exist between the samples and v^ in the absence 
of the noise voltages. As in section 3. 1. 2 we shall let be uniformly distributed, and 

a i i /2 

suppose that there is no fading replacing e by . 

In Appendix B (sec. B, 4) we show that the error in the differential phase estimate is 

A 8 

Variance <j> > . (163) 

n(4d|a 2 +4D2A 3 +D2A 4 ) 

The results of the last four discussions are summarized in Table 3, and the form of 
the bound is sketched in Fig. 21. 

As in the amplitude estimation case, the dependence upon the signal strength A is 
divided into four regions. When A is very small, the dependence is inverse with A, 
reflecting the fact that the receiver considers only the mean component in the samples. 
At higher values of A, we find that the performance is under -bounded by an inverse 
quadratic form. This form was obtained by assuming that a noisy phase reference signal 
was available. When we consider receiver structures, we shall find that such a sig- 
nal is obtained by using a second interferometer tuned to a slightly different spatial fre- 
quency, but occupying almost the same two aperture correlation areas. The error in 
estimating the phase difference between the two spatial frequencies then follows the form 
of this quadratic portion of the bound when we use a cross correlation type of receiver. 

As we shall see eventually, the bound in region 2 reflects the behavior of a receiver 
that tracks the phase of the fringes and in which the noise is so weak that it simply adds 
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linearly to the phase -tracking signal. Finally, in region 1 we have the performance of 
a receiver that extracts a perfect phase measurement so that the only disturbance that 
must be overcome is the random fluctuations in this phase caused by the atmosphere. 


3. 4 RECEIVER STRUCTURES FOR PHASE ESTIMATION 

Almost the same three receiver structures that were found to be optimum in the 
amplitude estimation problem apply also to phase estimation; that is, a fringe tracker 
when the signal strength is high, a- cross correlator for' intermediate signal levels, and 
an integrator when the signal strength is very weak. 

3.4.1 Fringe Tracker 

When the noise is small the instantaneous phase of the fringes can be accurately 
measured, and a phase estimate formed as 

N 

^ Im {In r . ) . (164) 

i=l 

When the noise is completely negligible, this estimate is, in fact, maximum- likelihood, 
and has a mean- square error 

2 

[$-<j>] 2 = -^-. (165) 


This is the same as the C-R limit of (158). As long as the noise is weak and a linear 
expansion of the logarithmic term is justified, the variance of an estimate based on 
noisy measurements is 



(166) 


where, as in the amplitude estimation problem, we have assumed that N is sufficiently 
large that the error is already small (10% or less). 

The effect of the nonlinear interactions of the noise are much more serious in the 
phase -estimation problem than in amplitude estimation. In order to exploit the zero- 
mean behavior of the phase disturbance, the phase cannot be measured in simply a 
modulo 2v representation; rather, the full value must be retained. If an instantaneous 
noise sample is larger than the signal component, the phasor may appear to encircle 
the origin, and an anomalous jump of 2 it rad may appear in the output phase -tracking 
voltage. This anomaly generates a gross error in the final phase estimate, as much as 
v rad if it occurs near the middle of the record. 

For this reason, the phase -tracking system is reliable only when the signal -to-noise 
ratio is above a threshold value. This is familiar in communication engineering as the 
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anomaly problem Inherent in nonlinear modulation; for example, the "click 11 problem 
in frequency and phase modulation. 


3, 4, 2 Integrator Receiver 

At the other extreme of very small signal strength, we again find that a receiver that 
considers only the mean component in the signal is optimum, and we form our estimate 
as 


A 

<t> - Phase 




< 167 ) 


When the ensemble average value of the r. and the measured average are fairly close 
(N large), the phase error is approximately 
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i N 

w 2 r i - r 

N i=l 1 


(168) 



The me an- square value then is 
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Again, under the assumption that is large, (169) becomes simply 
D„ 


[<m>] 2 *^- 


' z -+ 2 


b 2 d 2 ab 2 d 2 


When the signal is weak we obtain 


(169) 
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(171) 



±_ 

N 


? ? 1 

ab i d i 



This differs by only a factor of 2 from the 'C-R limit (159). Our approximation proba- 
bility overestimates the phase error by just this factor of 2, since only one of the orthog- 

N _ 

onal components of the difference (l/N) 2 r. - r. contributes to the phase error. Thus 

i=l 1 1 

it appears that for small values of A, the integrator is asymptotically efficient with 
increasing N. 

At high values of A {A > 2/D^), when the shot-noise becomes insignificant the error 
saturates as 





D, 


d i b i' 


(172) 


This is much larger than the saturation level for the phase -tracking system which 
achieves the C-R limit. This is due to the zero -mean component which remains in the 
signal and acts just as if it were a large noise voltage. 


3. 4. 3 Cross correlator Receiver 


Between the threshold values of A, below which the fringe tracker flounders, and 
the very small values of A where we know an integrator is optimum, is the last region 
that we must study. We know from our C-R limit that nothing better than l/A behavior 
will be possible. To even obtain this bound, we had to assume that an additional noisy 
phase reference signal was available and, with only the output of a single interferometer 
available, it would seem impossible to find anything better than the integrator receiver, 
once the signal level drops below the threshold value. 

It seems inconceivable, however, that we would ever want an estimate of the phase 
of only a single spatial Fourier component in the object. Certainly, if we attempt to 
reconstruct a full picture we need all of the Fourier components. It is then reasonable 
to suppose that we would have available the outputs of a number of interferometers tuned 
to different spatial frequencies. 

Suppose then that we look at two simultaneous sets of samples r^ and v^ of the form 


r - A 1 / 2 
i A 1 




a i + J b i 
e 1 + n. 




j«f> 2 cAjcL 


+ 1 *^, 


(173) 


where m and mu are independent shot-noise terms (independent because of orthogonality 
of different spatial frequencies). We shall assume that the atmospheric disturbances 
a i + jbj. and Cj + jd. are partially correlated. That is, the two interferometers share part 
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of the same aperture space. If we now crosscorrelate the two sets of data, we obtain 
a random variable R with an ensemble average value 
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R = 
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* A 

r i v i = A 1 


= A lA A l/2 
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(174) 


This average complex cross correlation has a phase that is simply the phase difference 
between the two spatial frequency components. An estimate of this phase difference 

Aj 2 = ^1 " ^2 can mac ^ e: 

N 

£ = Phase of — V r-v., (175) 

^ tvf Z_y 1 1 

i=l 

When N is large (error is small) we approximate the error in the same manner as in 
section 3, 4. 2 and obtain 


— fi — J5 (K - R ' Z 

: 1,2 LA 1,2 A 1 , 2 J ~ ( ^ } 2 • 


(1?6) 


This is the same technique that was used to estimate the error in the integrator 
receiver, and we find that 
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where 


* n(a.+c.) 

D n =e 


and 


(178) 


* jn(b -d ) 
B n = e 


If we regard AjA^ as the effective signal strength, we can see again that there are 
essentially three different behavioral domains. When the signal strength is very large, 
noise is negligible, and the phase error goes as 
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( 179 ) * 
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If the phases b. and are uncorrelated, their difference will be large and the 

denominator term will be extremely small. The phase estimate is then very poor, 
even when noise is absent, and it is of the same order as for the integrator receiver 
(172). Only the mean components of r. and v. are used by the receiver. 

If, however, we choose two neighboring spatial frequencies, the aperture displace- 
ments are almost the same, and we can use almost the same pair of aperture openings. 
The phase disturbances b- and d* will now be strongly correlated. Thus will be of 
the order of unity, and the phase variance will be much smaller, (See sketch below.) 



At somewhat lower values of effective signal strength the noise begins to influence 
the result, and the behavior is inversely with signal strength (assuming « A^), and 
we have 


r i 2 1 
^1,2^ = N 
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Finally, when the signal strength is very weak, we have an inverse quadratic 
behavior 
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When the phases b^ and are strongly correlated the result is similar to that predicted 
by the C-R inequality (163). 

Thus far, however, we have only considered the estimation of the phase difference 
between pairs of spatial frequencies. We can generate the absolute phase of the trans- 
form by an iterative procedure. Starting with the zero -frequency component (which must 
always have a zero phase), we make a set of differential phase estimates between pairs 
of neighboring spatial frequencies. These differential estimates are then summed to 
give the composite estimate. That is, 


Phase 0 v (k) = $ v (k) = A 0) 


+ A l, 2 + A 2, 3 


+ A- 


M-l, M* 


(182) 


where M is the total number of iterations. 

Everywhere but in region 1, the errors in these estimates are due primarily to the 
shot-noise terms, and hence all the differential errors will be uncorrelated. The 
total mean-square error is the sum of all of the mean-square differential errors given 
by Eq. 178. That is, 

M 

[$ v ®-+ v ®] 2 = X < 183 > 

m=l 

If we assume that all successive spatial frequency pairs are separated by the same 
amount, the amplitude and phase correlation terms will all be about the same. Further- 

i 

more, if the amplitude spectrum of the spatial transform is fairly uniform (all about 
the same), the successive errors will be of the same order, and the resultant phase 

error in © v (k) will go as M[e] , where M is proportional to the magnitude of k, and 
is given by 



(184) 


where Ad is the difference in displacement for neighboring interferometers. 

There will be some best choice for the interspatial frequency spacing Ad. When we 
choose a very small value of Ad, then the phase disturbances will be strongly corre- 
lated, but M will be large. Likewise, if we use a course spacing, M will be small, but 
the differential errors will be large. 

The mean-square errors for the three receiver forms, fringe tracker, crosscorre- 
lator, and integrator are summarized in Table 4 and the functional form is sketched in 
Fig. 22. The general conclusion is that when the noise is negligible, a fringe tracker 
is optimum and achieves the Cramer-Rao limit. At the other extreme, very weak signal, 


65 



Table 3. 


A 

Cramer-Rao lower bound to variance <j>. 
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Receiver characteristics for 




CRAMER - RAO LOWER BOUND 
Variance 0 ^ (l/N) G(A) 

REGION 1 
noise = 0 
strict 
lower bound 
for all A 

G X (A) = Cf b 2 

REGION 2 
a,b known 
strict 
lower bound 
for all A 

G 2 (A) = / 

ad 2 

REGION 3 
a known 
b uniform 
noisy 
reference 
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the integrator receiver is asymptotically efficient. In the intermediate region the cross - 
correlator receiver is best. Note, however, that increasing values of spatial frequency 
tend to shift this section of the performance curve upward, because of the increased 
number of iterations that are required* This tends to increase the value of A that 
determines the choice of the integrator or crosscorrelator receivers. Indeed, if kk is 
large enough, it would appear that the cross correlator is no longer useful, and we use 
either the fringe tracker or integrator, depending upon the strength of the signal -to- 
noise ratio relative to the threshold value. 

We should also note that the iteration technique may flounder when there is a gap in 
the spatial transform of the object. 

3.5 INTERPRETATION OF RESULTS 

3.5.1 Extension to Samples in Other Variables 

We have dealt entirely with time samples of the mutual-coherence process. We now 
want to extend our results to the general case in which we have samples in other 
domains. When we have additional samples for aperture position, these samples are 
essentially independent, and should be treated in the same manner as the time samples. 
That is, we should simply integrate across both time and space when the signal level 
is very low, take square magnitudes and then integrate when the signal level is in the 
intermediate region, and average the phase and log amplitude when the signal-to -noise 
ratio is very high. 

If we know that the object has essentially the same pattern over a broad range of 
wavelengths, or we want the color -averaged pattern, we treat the optical-frequency 
samples in the same manner. 

If a full picture of the object is required, then we must estimate each spatial- 
frequency component from the corresponding displacement samples. The interdepen- 
dence between samples of different spatial frequency can be exploited to some extent. 

For example, our iterative phase -estimation procedure uses such statistical dependen- 
cies. We also suspect that a similar crosscorrelation procedure would aid in estimating 
the amplitude of a weak spatial-frequency component if we have nearby a strong compo- 
nent that in effect could be used as a channel phase -reference signal. 

3.5.2 Numerical Examples 

Thus far, our results are still couched in the parameters A, B^, and D^. In order 

to assess the applicability of the estimating structures that we have obtained, we turn 

now to consider practical values of these parameters. Let us consider Fig. 20, which 

describes the mean-square error for various receiver forms in estimating the magnitude 

of the spatial transform. As an example, suppose the variances of the log-amplitude 

2 2 2 

fading and the phase processes are ^ 0.5 and » 50 rad . These values were found 
typical for a horizontal 4. 5 km path (see Section IV). With these values, the curve of 
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-210 dB 10 log 1Q A -7 dB +7 dB 


Fig. 23. Receiver characteristics. 

Fig. 20 takes the form shown in Fig. 23. 

The most immediate reaction to Fig. 23 is the extremely large range over which 
region 3 seems to apply. Over more than a 200-dB spread, it appears that square -law 
(or cross correlator) receivers would be best. 

In order to appreciate this figure, we must have some idea of the order of magnitude 
of signal-to -noise ratio to be expected. Consider, for example, the Sun, which has an 
approximate intensity level of 

1. 8 X lO" 1 W/(cm) 2 sterad A 0 . 

For typical atmospheric parameters, let us suppose 

Ax 2 = 1 (cm) 2 (10~ 4 (m) 2 ) 

AT - 1 msec (10 -3 sec) 

Av - 50 A". 

o 

We would then have incident on our detector approximately 9 X 10 j/sterad. 

-19 

Since the energy per photon is approximately 10 J, this would be detected 

as 9 X 10 -1 ^ photons/sterad, and within just the diffraction-limited disk for a 1-cm 

—12 8 
aperture (0. 25 X 10 sterad), the resultant photon count would be ~2 X 10 photons. 

7 

For 10% quantum efficiency this becomes 10 photoelectrons. If the mean-square 

value of the fringes is approximately 10% of the average squared signal level, we end 

6 

up with a signal -to-noise ratio of the order of 10 , or 60 dB. 

From Fig. 23 we see that this value falls in region 1, in which a fringe -tracking 
receiver would be best. For objects that reflect sunlight, such as the Moon, the 
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intensity level drops by a factor of ~10 (70 dB). Under the same conditions as before, 

this puts our signal -to -noise ratio down to approximately -10 dB, well within the region 
for the square-law receiver. 

3. 5. 3 Comparison with Existing Techniques 

We shall now discuss briefly the comparison of the preceding receiver structures 
with two existing techniques. 

a. Long-Exposure Photography 

The first of these is the technique of long-exposure photography followed by inverse 
filtering. This long-exposure technique is, in essence, an integration of the mutual- 
coherence process in both time and space over many correlation regions. Thus it is a 
realization of the receiver structure whose performance is denoted by in Figs. 20 
and 22. As we can see from our recent results, for many objects this technique is quite 
inferior to the square -law and tracking receivers. It does possess many practical 
advantages, however, probably the most important being the ease of implementation. 

All of the integration can be done on a single plate, and many spatial frequencies can 
be detected simultaneously. 

b. Hanbury Brown and Twiss Intensity Interferometry 

Another technique for estimating the amplitude of the spatial-frequency components, 

11 12 

is by correlating the intensity of the field at two points. ' For Gaussian processes, 
the fourth- and second-order moments are conveniently related. Thus, from knowledge of 
the fourth-order moment, the second-order or mutual- coherence function can be deter- 
mined. The main disadvantage with this technique is that large detector bandwidths are 
required if we are to recover any significant amount of the intensity fluctuation in the 
light. This detector bandwidth turns out to be the limiting factor, the bandwidth of the 
light itself being of less significance. It is interesting to compare the performance of 

this scheme with that for the square -law processor. Given that both schemes operate 

2 

within the correlation bandwidth Av, with aperture areas Ax , the ratio of mean-square 
errors turns out to be roughly 


Id 


2 

square law 


Id 


2 

Brown and Twiss 


BW 
7 9 

(Av) Z 


where B is the bandwidth of the two detectors employed in the intensity correlation, and 

W is the bandwidth of the atmospheric fluctuations. Generally the product BW is much 
2 

less than (Av) , and the Brown and Twiss method would be inferior. We assumed, how- 

2 

ever, that the receiving areas were matched to the correlation area Ax . If the entire 
available aperture area is A, the signal -to -noise ratio in each detector output would be 
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increased by the factor A/Ax , with the result that the mean-square error for the Brown 
and Twiss method would be reduced by the factor [Ax' /A] . On the other hand, with the 
square-law receiver, simply increasing the area does nothing to the mean-square error 
unless we employ parallel processing of the aperture samples. In this case, the error 
decreases by the factor Ax /A, which is the reciprocal of the number of essentially 
independent spatial samples. Thus the Brown and Twiss system is able to make more 

2 , 

efficient use of the aperture. Nevertheless, because of the great disparity between (Av) 
and BW, the system still is inferior to square-law processing of the mutual-coherence 
samples unless the aperture is very large. 
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IV. CHANNEL MEASUREMENTS 


An experimental investigation was conducted in order to measure some of the prop- 
erties of the phase and amplitude disturbances induced by the turbulent atmosphere which 
play a part in the high -re solution imaging problem. The purpose of these measurements 
was not only to obtain channel data but, perhaps more important, to gain practical 
insight,- appreciation, and respect for the turbulent atmosphere as a channel and for the 
problems of designing equipment for use on this channel. We now have the twofold pur- 
pose of (i) illustrating the nature of the atmospheric channel, and (ii) describing and dis- 
cussing some of the ways in which the channel can be probed. 

The optical path employed for most of this study was between The Museum of 
Science, in Boston, and the Harvard College Observatory, in Cambridge. The map in 
Fig. 24 shows the locations of these sites which are separated approximately 4. 5 km. 

The light source was a Helium-Neon laser housed in a weathertight box located on the 
tower balcony of the museum as shown in Fig. 25. The transmitted light was collimated 
into a beam, 2 inches in diameter, having uniform phase and approximately Gaussian 
intensity profile. A small spotting telescope mounted at the rear of the laser was used 
in the initial alignment of the beam with the receiver. Precise pointing of the beam 
required persons at each end of the link equipped with radio telephones for communica- 
tion. 

Fastened to the top of the box there is a smaller case containing a mercury arc lamp 
of the kind used in street lighting. By placing various gratings in front of the lamp some 
of the qualitative effects of the turbulence upon image degradation were studied visually. 

The roof of the Harvard College Observatory was the terminal point of our link. We 
were privileged to be able to make many of the visual and photographic observations with 
the famous 15-inch refracting telescope shown in Fig. 26. Although it is more than 100 
years old, this telescope is still in excellent condition, and it is only the growth of the 
metropolitan area with subsequent deterioration of observing conditions which has 
reduced its usefulness for astronomical observing. This telescope is ideal for our pur- 
poses, as it has a large unobstructed aperture and long focal length. The advantage of 
long focal length is high magnification in the primary focal plane. 

Because some of our experimental electronic measuring equipment was too bulky and 
awkward to use with the 15 -inch telescope, an 8-inch reflector was mounted on a bench 
in an adjacent penthouse. 

4. 1 QUALITATIVE CHANNEL DESCRIPTION 

4. 1. 1 Aperture Intensity Profile 

The property of the disturbance that is most easily measured is the random intensity 
profile of the incident wave. In a plane that is a few centimeters ahead of the focal plane 
of the telescope, the field, although spherically converging, still retains the same 
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Fig. 24. Optical path. 



Fig. 25. Laser (0. 6328 \x) atop The Museum of Science, Boston, Massachusetts. 



Fig. 26. 15 -inch refractor at Harvard College Observatory, 

C ambridge , Massachusetts . 
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intensity profile that is incident on the objective. The linear scale is conveniently 
reduced, however. Either by placing the photographic film in this plane or using an 
additional lens to project this plane on the film, the spatial variations in the incident 
intensity are recorded. Figure 27 shows 8 successive frames taken in this manner with 
a 16-mm motion picture camera coupled to the 15-inch refractor. The aperture corre- 
lation region for the amplitude disturbance can be roughly estimated. In the example 



Fig. 27. Intensity profile across the 15 -inch objective. 


shown here it appears that intensity correlation is lost beyond approximately 1 inch. 
Some of the time behavior can also be studied with this motion-picture technique. When 
the film is run on a motion projector, the effect of the wind blowing the disturbances 
across the aperture is very marked. The hypothesis that the disturbances are "frozen" 
into the atmosphere does not hold indefinitely, and from the film it appeared that after 
3 inches or so of movement across the aperture the structure had greatly decayed. Such 
joint time and space correlations of the wavefront intensity have been studied by J. E. 

l ft 

Ehrenberg and the results have been published in his Master's thesis. 

4. 1. 2 Focal-Plane Patterns for Small Apertures 

The patterns of Fig. 27 describe only the intensity variations of the field. For both 
digital communication and image -estimation systems, however, the phase of the process 
is usually of greater interest. Some of this phase information was deduced simply by 
varying the aperture size and observing the behavior of the focal-plane image. When 
we employ a circular aperture that is smaller than both the phase and amplitude corre- 
lation scales, the image in the focal plane is an almost perfect Airy pattern. That is, 
it is the true diffraction pattern from a circular aperture, showing a central disk with 
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Fig. 28. 

Focussed images, 0. 5 -inch 
diameter pupils. 


Fig. 29. 

Focussed images, 1-inch 
diameter pupils. 


Individual 

Diffraction 

Patterns 



Complete 

Diffraction 

Pattern 



Fig. 30. Aperture sampling of the image. 




a series of increasingly fainter concentric rings. The only effect of the turbulence is 
then a time -variant fluctuation in the over-all intensity of the pattern. Variations in the 
absolute phase of the wave are undetectable in such a recording, but are of no interest 
in the imaging problem. As the diameter of the aperture is slowly increased, the shape 
of the pattern remains, but the dimensions shrink because of the increased aperture 
size. The first noticeable effect as the correlation diameter is exceeded is a random 
motion of the Airy disk. This random motion can be related to the first-order or linear 
term in a Taylor* s series expansion for the phase disturbance, and provides an almost 
complete description of the process when the aperture is small. Finally, when the aper- 
ture is made much larger than the correlation scale, the simple structure of the distur- 
bance is lost, and a highly distorted image is formed. 

Some of these effects can be seen in Figs. 28 and 29. These pictures were recorded 
on the same evening as the previous intensity measurements and were obtained in the 
following manner. A mask consisting of an array of small holes was placed over the 
objective. The focal-plane pattern, which comprises the diffraction pattern of this entire 
array, is not too informative. By observing the pattern just a bit forward of the true 
focal plane, however, we can resolve the focussed image of each aperture hole sepa- 
rately, as illustrated in Fig. 30. 

Two sequences of motion- picture film were shot in this manner: First with 0. 5-inch 
diameter holes on a 2 X 2 inch grid, and then with 1-inch diameter holes on a 4 X 4 inch 
grid. The first example (0. 5-inch diameter holes) gave Airy disk patterns showing 
almost no discernible motion when the film strip was viewed with a projector. It would 
appear that over the area of these 0. 5 -inch diameter holes, the field remained quite 
coherent in both phase and amplitude. The second set of pictures, made with 1-inch 
diameter holes, showed appreciable motion, perhaps half the disk width which corre- 
sponds to a wavefront tilting of 15-20 prad. These two sets of data illustrate the point 

1 9 

brought out by Fried, that random tilt is a very significant characteristic of the phase 
disturbance over a small aperture and that if this tilting or image motion can be tracked, 
the useful aperture can perhaps be increased approximately three times before the 
higher order disturbances significantly distort the image. 

It is obvious, however, that even in the 1-inch diameter case, the tilting does not 

completely characterize the disturbance. Already there is some loss in the symmetry 

of the pattern, caused by either spatial variation in the intensity or the higher order 

phase disturbances. J. H. Shapiro has carried out an investigation of this problem, 

seeking simple models to characterize the disturbances over small apertures. His 

20 

results are published in his Master's thesis. 

4. 1. 3 Focal-Plane Patterns for Large Apertures 

When the full aperture of the telescope is employed, the focal-plane pattern no longer 
has a simple structure. Figure 31 is an example of the pattern resulting when the full 
aperture of the 8 -inch reflector was used, and about all that can be said is that the 
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Fig. 31. Focussed image for full 8-inch objective. 




Fig. 32. Profile of intensity across 8-inch objective. 
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energy in the image is randomly distributed over approximately 40 ^rad field of view. 

The incident intensity profile over the aperture was also measured a few moments 
before, and a sequence of frames is shown in Fig. 32. The similarity of this intensity 
pattern with that of Fig. 27 suggests that conditions were somewhat the same on the two 
different nights. The dark disk in the center of each pattern and the rings on each side 
are simply the shadows of the secondary mirror and supporting structure. 

4. 1. 4 Fresnel Biprism Measurements of Field Coherence 

More information concerning the phase process can be obtained by using the Fresnel 
biprism form of interferometer described in Section II. This was set up with a motion 
picture camera as shown in Fig. 33. 



Film 


Fringes are formed in the region where the two displaced halves of the aperture 
overlap. Recall that the spatial phase and intensity of these fringes is directly related 
to the phase and magnitude of the incident mutual coherence function. The distance 
between the prism and film determines the amount of overlap of the aperture sections, 
and hence the displacement variable of the coherence measurement. Examples of this 
technique are shown in Figs. 34 and 35, in which two different displacements of 2 inches 
and 4 inches were used. The fringes of interest run horizontally, and are the most 
closely spaced ones that are visible. The other coarser fringes result from diffraction 
from the boundaries of the aperture and the discontinuity in the prism. 

If a uniform plane wave had been incident, the interference lines would have been 
absolutely straight. The presence of corrugations in the incident wavefront causes the 
fringes to wiggle, and a particular fringe line follows the locus of constant phase dif- 
ference between points of given displacement. 

Various forms of distortion can be deduced by the character of the fringe pattern. 
Uniform tilting of the wavefront across the entire aperture is undetectable, but can be 
discerned with other methods (see section 4. 1.2). A differing tilt on the two areas that 
are interfered is observed as a local twisting, bunching or stretching of the fringes. 
One can obtain a spatially continuous measurement of the distortion by following along 
the fringes, whereas in the Airy disk method the disturbance is only coarsely sampled. 
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Admittedly, the distortion in the fringes is hard to see in these small prints. When 
the film is projected with high magnification onto a screen, however, twisting of the 
fringes is quite noticeable in both examples. When we compare these patterns with the 
interferogram for an undistorted incident wave, we find that the two sets of lines are 



Fig. 34. Fresnel biprism pattern — 
2 -inch overlap. 



Fig. 35. Fresnel biprism pattern — 
4-inch overlap. 


correlated over distances of only 1 inch or so. The significance of this to image esti- 
mation is that there are areas whose fringes differ in phase by approximately tt radians. 
These will tend to give a cancellation effect if the mutual-coherence function is inte- 
grated over an area larger than the correlation region. The average signal -to-shot-noise 
level will not improve if the integration area is increased beyond this value, but we lose 
some of the channel diversity that could be exploited if proper spatial sampling and one 
of the post processing schemes discussed in Section 111 were used. 

The effect of the atmospheric fading on the mutual coherence is also brought out in 
these two pictures. We recall that if the amplitude on either of the points being inter- 
fered should fade, then the coherence function would fade, since it is the product of the 
fields. Thus the strongly fringed regions appear as clumps in the interferograms. 
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4. 1. 5 Full-Aperture Phase Mapping — Spatial Heterodyning 

This Fresnel biprism technique yields relative measurements of the phase distur- 
bance between points. Often it is desirable to have a full map of the absolute phase pro- 
cess. Probably the nicest picture that can be made of the phase across the aperture is 
obtained by beating the incident wave against a tilted, but otherwise uniform, plane wave 
of the same frequency. At optical wavelengths this reference wave must usually be 



Fig. 36. Spatial heterodyne pattern. 


obtained from the same source as the incident wave in order to retain temporal coher- 
21 

ence. M. Tamny has studied such a measurement technique, using a two-way path, 
reflecting the wave from a downrange mirror. Figure 36 is one of the pictures described 
in his Master's thesis, and shows the type of pattern to be expected. This was obtained 
over a relatively short path of a few hundred feet for which the amplitude disturbances 
are minimal. The phase disturbance across the aperture shows up in the twisting of the 
fringes. 

4.2 QUANTITATIVE MEASUREMENTS OF PHASE AND AMPLITUDE 
- TWO-POINT COHERENCE MEASUREMENTS 

This, then, is the extent of the photographic measurement techniques that were 
studied. They yield reasonably qualitative pictures of the turbulence -induced distur- 
bances. The most important property of these disturbances with respect to image esti- 
mation is often the temporal variations in the field coherence between two points. 

4.2.1 Description of Equipment 

In order to obtain a more quantitative understanding of these properties, the equip- 
ment pictured in Fig. 37 was employed. Here we used the telescope as an interferom- 
eter, as described in Section H. The aperture was completely masked, except for two 
holes at the points where we wanted to measure the field coherence. The reach of the 
interferometer was extended by means of the mirrors (M) in Fig. 38. The diffraction 
pattern for this aperture mask is an Airy disk crossed by straight fringes, whose 
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Fig. 37. 8-inch reflector with aperture mask and fringe scanner. 



Fig. 38. Extension using mirrors. 
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spatial phase and amplitude correspond to the phase and magnitude of the field coher- 
ence function. 

a. Fringe Scanner 

These fringes were scanned with the device illustrated schematically in Fig. 39a and 
pictured in Fig. 39b and 39c. 

The focal plane -field pattern was projected by lens L j onto the Ronchi ruling held 
by the sliding block F. This block was pushed back and forth by means of the motor, 
cam assembly C, and a return spring. As the Ronchi ruling moved back and forth the 
fringes were scanned by the lines in the ruling, and the transmitted light energy was 
focussed by lens onto the photodiode detector D^. The signal from this detector con- 
sisted of: (i) a pair of low-frequency terms corresponding to the fluctuations in the total 
incident energy through each hole, and (ii) a periodic signal with phase and amplitude 
related to the phase and amplitude of the fringes. With the particular cam shown in the 
picture, the Ronchi ruling moved back and forth twice each second, a total travel of 
2 inches per second. Since the Ronchi ruling contained 300 pairs of alternate opaque 
and transparent lines per inch, the mean frequency of the periodic output signal was 
600 Hz. 

Another cam, which moved the ruling at only half this speed, was made, and it 
reduced the carrier frequency to 300 Hz. It is always necessary, however, to use a 
carrier of high enough frequency to avoid aliasing the interference signal with high- 
frequency components in the baseband intensity fading signals. On the other hand, the 
higher frequency cams have the disadvantage that there is a shorter period between turn 
around moments for the ruling. When the ruling stops and starts to go in the other 
direction, it is usually impossible to retain the phase of the signal; thus, the useful 
observation interval is limited. 

The actual phase of the fringes in the diffraction pattern was obtained by comparing 
the phase of the signal from D ^ with that of a reference signal generated by a small 
lamp A, lenses and L^, and detector D^. The lamp contained a small helically 
wound filament that was focussed onto another section of the Ronchi ruling, with the 
spacings on the helix matched to the lines on the ruling. 

When the phases of these two signals were compared, the effect of variations in the 
velocity of the ruling cancelled out, and we obtained the phase disturbances that were 
due only to the atmosphere. 

The block diagrams of the recording and processing instrumentation are shown in 
Fig. 40. Both of the signals were filtered to remove low-frequency components, and 
then recorded on a dual-track magnetic tape recroder. Comparison of the two signals 
was made later in the laboratory, by using a pair of demodulators to obtain the in-phase 
and quadrature components of the wave coherence function. These two voltages, together 
with the signal envelope, were then sampled, converted to digital form, and re-recorded 
for playback later into the TX-0 digital computer. With the computer we then calculated 
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Data Recording 



the phase angle of the coherence process, and derived various statistical properties of 
the phase and amplitude disturbances in the wave coherence function. 

b. Alignment 

The procedure for adjustment of the scanning device was the following. The required 
projection distance between the lens and the Ronchi ruling was calculated, and final 
focussing adjustments were made by observing the output from and adjusting to obtain 
maximum amplitude in the periodic component. The scanning device was then tilted, or 
the aperture mask rotated, to obtain parallelism between the fringes and the lines on the 
ruling, again adjusting for a maximum signal amplitude. The photodetectors D 1 and 
were mounted on plates that could be moved about, so that the light beams that exit from 

and were centered on the small detecting surface of each diode. 
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Fig. 42. Signal waveform 
(25 msec/cm). 


Fig. 43. Signal waveform 
(5 msec/cm). 
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Use of mirrors to extend the reach of the interferometer brings up further alignment 
problems. The telescope tube was first rotated so that one could look into the eyepiece 
(lens L^). Without the center mirror in place, this eyepiece focussed to give a proper 
image of the source to the eye (virtual images at infinity with the two Airy disks 
appearing to completely overlap). With the center mirror installed and the mirror sup- 
ports adjusted, the two Airy disks again overlapped in the same place. This established 
a proper relationship between the optical axis of each pair of mirrors and the telescope. 
The procedure for alignment of the scanning device was then as before. 

c. Typical Output Signals 

An example of the fringe pattern that was produced in the focal plane is shown in 
Fig. 41. The hole diameter should be so chosen that the field remains coherent across 
the area of the hole. That is, so that there is little noticeable movement of the Airy 
disk. In the example shown, we have intentionally made the hole diameter a bit larger 
than is advisable, in order to illustrate an important point. The Airy disks tend to sep- 
arate if the holes are too large, because of differences in the local tilt or angle of 
arrival. The common area in which the fringes are generated will then vary randomly. 

In such case, fading in the mutual -coherence process results not only from the ampli- 
tude fading of the waves but also from the phase incoherence. 

Figures 42 and 43 are typical of the waveforms recorded from the signal and refer- 
ence detectors. Figure 42 has the longer horizontal sweep scale, and shows the inter- 
ruptions in the reference signal occurring at the instants when the direction of the ftonchi 
ruling is reversed. These interruptions provided a convenient signal for triggering syn- 
chronization circuits in the processing equipment, Figure 43 shows a portion of the 
same record on an expanded scale. The fading that often occurs in the interference 
signal is well illustrated. 

A typical plot of the phase of the process as calculated by the TX-0 computer is 
shown in Fig. 44. 

4. 2. 2 Springtime Results — 4. 5 -km Path 

The first results that we present were obtained on the Museum of Science -Harvard 
path, late in April 1968. Experiments were conducted on several evenings and the par- 
ticular results presented here were representative of this set of experiments. 

These particular recordings were taken on April 28th, 1968 between 1 a. m. and 
3 a. m. under the weather conditions tabulated below. 


air pressure 

1017,0 mbar 

temperature 

48*F 

dew point 

38 “F 

wind bearing 

330° 

wind velocity 

1 2 knots 

velocity perpendicular to path 

4 m/sec 
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Hole diameters of 0. 5 inch were necessary in order to get within the correlation 
scale of the disturbances. 

Recordings of the time -variant mutual- coherence process were made for a range 
of displacements, and the resultant data were processed for (i) the probability distribu- 
tion function of the amplitude of the mutual- coherence process, and (ii) the probability 
density function and temporal and spatial structure functions of the phase of the mutual- 
coherence process. 

a. Amplitude Statistics 

The distribution function of the amplitude process is shown in Fig. 45 for displace- 
ments of 2, 4, and 6 inches. The results are plotted with log-normal probability coor- 
dinates. A log-normal distribution function then appears as a straight line on such a 
scale system. The three curves have been arbitrarily scaled in amplitude to separate 
them on the graph. 



0.1 99.8 

Fig. 45. Probability amplitude distribution function. 

The first point to note is the fit of the data to straight lines which substantiates the 
hypothesis that the mutual-coherence process suffers log-normal fading in propagating 
through the atmosphere. 

The severity of log-normal fading is described by the variance of the log amplitude. 
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This is the slope of the distribution function when plotted with the special coordinates. 
The standard deviation, a, of the natural (rather than base 10) log amplitude of the 
coherence process is tabulated below . 


Displacement 

Standard Deviation 

2 Inches 

cr “ 0 . 67 

4 inches 

Q 

n 

o 

H-* 

6 inches 

ff = o. 64 


b. Phase Statistics 

Measurements of the phase process were conducted in the following manner. Each 
half-cycle of the Ronchi ruling was taken to constitute one record of the phase process. 
Depending upon which cam was used (300 Hz or 600 Hz) this record had a useful length 
of 500 or 250 msec. At the end of the record, the Ronchi ruling stops and reverses its 
motion, which makes it rather difficult to retain the continuous value of the phase across 
this interruption. Indeed, the change in direction of the ruling reverses the sign of the 
phase between the fringe and reference signals. 



0.1 99.8 


Fig. 46. Probability phase distribution function. 
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In processing each record, the computer arbitrarily set the phase of the first sample 
to a zero 2ir modularity. From then until the end of the record the full (rather than 2ir 
modulo) value of the phase was retained. The computer next calculated a set of phase 
differences between the initial sample and samples of increasing delay in the record. 
From a large number of records, the probability amplitude distribution function was 
then determined. A typical set of such data is plotted in Fig, 46. The scales were 
selected to yield a straight line for a Gaussian distribution function, and again the 
Gaussian hypothesis appears reasonable. 

From the slopes of these curves we obtain the variance of the phase difference, and 
hence the temporal structure function of the phase of the mutual- coherence process. 
These functions appear in Fig. 47. We note that even with 4-msec delay the standard 
deviation in the phase difference is approximately 1 radian. By the end of the record, 
the variance of the phase difference appears to be saturating at approximately 12 radians. 
Assuming that when the structure function saturates the local variations in the phase 





Fig. 47. Temporal phase -structure function. 
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process are uncor related, we can obtain the variance of the difference in phase between 
the two aperture points simply by dividing the saturated value by two. This method of 
estimating the phase variance has the advantage that the mean value of the phase, which 
is a function of the relative positions of the lamp filament and diffractions images, need 
not be considered. All three curves appear to saturate at about this same value of 
12 rad, which implies that the standard deviation of the phase disturbance in the mutual- 
coherence process is IZ/\J~Z = 8.5 rad for all three displacements. 

It is interesting to note, in both the amplitude and phase measurements, an almost 
complete lack of dependence upon the displacement variable. Additional measurements 
that were made out to 18-inch displacements also showed this effect. This suggests that 
the spatial structure function of the amplitude and phase processes has saturated at dis- 
placements beyond even a very few inches, and hence that there is little correlation in 
the local wavefront disturbances. 

As we have mentioned, the data recorded in this series of late April experiments all 
showed about the same behavior, as the weather conditions remained quite similar. 
Attempts to obtain phase measurements earlier in the year were not entirely successful, 
because of the extremely rapid fading which tended to alias the 600-Hz carrier of the 
coherence process. 

4. 2. 3 Autumn Results — 4. 5 -km Path 

We do have one other set of measurements that were obtained on this path, and show 
the channel in a much more temperate mood. These measurements were obtained on 
the evening of October 25th, 1967 under the following weather conditions. 


air pressure 

1012. 1 mbar 

temperature 

57 °F 

dew point 

55°F 

wind bearing 

130° 

wind velocity 

15 knots 

velocity perpendicular to path 

2. 6 m/sec 


The conditions on this evening were such that a 1-inch diameter aperture was within 
the correlation scale of the turbulence, a good indication that the atmosphere was more 
stable than in the springtime experiments. 

Figure 48 shows the amplitude distribution function from which we find that a « 0.43, 
and is a condition of considerably less severe fading than in the former examples. 

The temporal phase -structure function is shown in Fig. 49. Here we find that the 
correlation time of the process is considerably longer than in the other measurements. 
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Furthermore the saturation value of 5 rad is only half that encountered in springtime. 
These autumn measurements were made under conditions of somewhat less 




Fig. 48. Probability amplitude distribution 
function. 


Fig. 49. Temporal phase -structure 
function. 


crosswind and much higher humidity than the spring measurements. Although this may 
be the reason for some of the difference, the longer term climate difference may have 
had an influence. 

4. 2. 4 100-m Roof Path 

Some phase measurements were also conducted over a much shorter path (100 m) 
along the roof of M. I. T. Building 20, in order to look at some of the differences in the 
disturbances between long and short paths. These measurements were made on 
April 14th, 1968 under the following weather conditions. 


air pressure 

1018. 5 mbar 

temperature 

45°F 

dew point 

27 Q F 

wind bearing 

070° 

wind velocity 

10 knots 

velocity perpendicular to path 

2. 2 m/sec 


Apertures of 0. 5-inch diameter were selected to get within the correlation scale. The 
amplitude of the coherence signal showed little significant fading, because of the short- 
ness of the path. The temporal phase -structure functions for these measurements are 
shown in Fig. 50. The notable feature of these measurements is the magnitudes of the 
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Fig. 50. Temporal phase -structure function. 

phase variances which are almost the same as those presented in section 4. 2. 2, although 
the path length is much less than that between The Museum of Science and Harvard Col- 
lege Observatory. This seems to support the contention that the phase disturbances are 
generated by the nearby turbulence, while the amplitude disturbances are a phenomenon 
of the long atmospheric paths. 

4. 2. 5 Laboratory Simulation 

A final set of measurements was conducted in order to evaluate laboratory simulation 
of a turbulent atmospheric path (Fig. 51). The fringe scanner and a small telescope were 
mounted on one bench and the laser and collimator on another. In front of the telescope 
objective we placed a strong source of local heating, an electric hot plate. A fan was 
located 5 feet away from the hot plate (directed somewhat away from the objective) in 
order to stir up the air in the room. The actual locations of the hot plate and fan, 
and the position of the equipment in the room all tended to influence the severity 
of the turbulence. By adjustment of the various components, the amplitude of the 
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Laser 


Telescope 





Fig. 52. Temporal phase -structure function. 


mutual-coherence function can be made to undergo fading that is as severe as that 
encountered on the long atmospheric paths. In one particular example phase measure- 
ments were made using holes l/8 inch in diameter. A set of temporal phase -structure 
functions is shown in Fig. 52. This shows much the same character as was obtained on 
paths through the atmosphere. The maximum variance in the phase is about the same 
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as obtained on the long path (springtime), although the correlation time seems to have 
about twice the former value* This last point can easily be taken into account in any 
simulation study* 

Both of these laboratory measurements and the measurements made on the roof show 
only a small increase in the phase variance as the displacement is extended. This would 
all seem to suggest that the spatial structure function of the phase of the wavefront tends 
to saturate at just a very few correlation diameters. (By correlation diameter we 
mean the maximum aperture diameter for which a well-formed stationary Airy disk 
is observed.) 

4. 3 DISCUSSION OF RESULTS 

As we have stated, our intention in this phase of the research was not to conduct an 
investigation in atmospheric physics but rather more to understand the general features 
of the turbulent atmosphere as a channel. The results obtained over the Museum of 
Science -Harvard College Observatory link are probably reasonably representative of 
the disturbances in the mutual-coherence function that are to be encountered over long 
horizontal paths, for example, as would be encountered in observing objects low on the 
horizon. 

The results for the roof-top path are probably somewhat more representative of the 
shorter paths near zenith, although the channel may be a somewhat severe example, 
since it is a roof top with several ventilators and other sources of agitation. 

We also demonstrated that a simple laboratory model can be used to give a reason- 
able simulation of the major effects of the turbulent atmosphere in distorting the 
received mutual- coherence process. 

Finally, we would like to make a few comments on the coherence -measuring tech- 
nique. It seems that the basic principle of scanning fringes is a good way to obtain the 
phase and amplitude processes. The actual scanning mechanism could be improved. The 
main problem that we encountered was inability to separate the interference signal from 
the baseband fading components when the fading was very rapid. The solutionis to use 
a carrier of, perhaps, 2-3 kHz. With available Ronchi ruling line densities the recip- 
rocating scanning techniques do not appear to be the best approach, however. A better 
approach may be to use a Ronchi ruling built into the rim of a scroll case which could 
be rotated at quite high speed. Another, and perhaps easier, approach would be to use 
a rotating disk consisting of lines ruled like the spokes on a wheel. If the diffraction 
pattern were confined to a small area near the periphery of the disk, the spokes would 
appear almost parallel, and the system would work just the same as for a conventional 
Ronchi ruling. A device of this latter form has been developed by Ramsey and Kobler*^ 
for monitoring the atmospheric seeing conditions. In all of these mechanical scanning 
techniques it appears necessary to maintain a reference signal such as the light-bulb 
arrangement used on our own scanner. 
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V. CONCLUSION 


5. 1 SUMMARY OF RESEARCH 

The object of this work was to develop a reasonably thorough understanding of the 
limitations and possibilities of image estimation through the turbulent atmosphere. 
Although some approximations were unavoidably introduced to simplify the modelling 
procedure, we believe that the results still reflect the salient features of the general 
case. 

We showed that if a receiver were given complete freedom in making field measure- 
ments, it could extract a sufficient set of data samples by determining the local mutual 
coherence of the incident radiation. That is, we found that the local mutual coherence 
of the field, rather than the instantaneous field function, constituted a statistically suf- 
ficient description of the process. Rather than extract a continuous record, a receiver 
may sample the process appropriately in the variables of frequency, time, position, and 
displacement (spatial frequency). We then studied the ways in which these coherence 
measurements are obtained and considered the unavoidable disturbances introduced by 
the detection device. 

The resultant samples were then shown to be well modelled as random variables 
having Gaussian log-amplitude and phase caused by the atmosphere, and additive Gauss- 
ian noise components injected by the photodetection process. The extraction from these 
noisy coherence measurements of estimates of phase and amplitude of the spatial Fou- 
rier components in the object was next considered. Lower bounds to the variance of 
parameter estimates that could be expected from any unbiased processor were then found 
to fall into four general regions. When the signal strength is very high and the detec- 
tor noise negligible, the variance is dependent only upon the severity of the turbulence. 

In this region a receiver that tracks the phase and log- amplitude distortion can achieve 
the lower bound. At lower signal levels the noise comes into play, and in region 2 it 
has a linear effect on the variance of both phase and amplitude estimates. This region 
is usually rather small, and lowering the signal-to-noise ratio brings us to region 3. 
Here the noise power acts quadratically on the variance of the amplitude estimate. A 
receiver that averages the square magnitude of the samples is almost optimum. The 
phase estimation problem, however, exhibits a threshold effect in this region. Once the 
noise is dominant in each sample there is no way to obtain good phase estimates with a 
sequence of samples at only one spatial frequency. If simultaneous measurements are 
available at neighboring spatial frequencies, estimates of the differential phase can be 
obtained whose variance has the quadratic form similar to that for amplitude estimation. 
Finally, when the signal-to-noise ratio is very poor, a fourth region may exist. Pro- 
vided the object and receiver are rigidly secured, the variance may again increase only 
linearly with the noise power, but the only part of the signal that is really of use is the 
small residual mean component. The best receiver in this case is one that simply 
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integrates the coherence measurements. 

The optimal techniques were then compared with various existing systems and some 
practical examples were considered. 

The results of experimental channel studies were presented. Photographic measure- 
ments prove quite illustrative of the turbulent effects on the wave -front coherence. Spa- 
tial and temporal structure functions were also determined on both long and short optical 
paths and provide a more quantitative understanding of the severity of the disturbances 
encountered. 

5. 2 DESIGN PHILOSOPHY 

From the examples presented, we are convinced that there is considerable merit in 
employing the more sophisticated data processing over the existing techniques. Indeed, 
it appears that a large class of problems will fall within regions 2 and 3, in which the 
square law and crosscorrelator receivers have truly significant gain over the integra- 
tion technique. 

Unfortunately, just as in channel coding and decoding, our enthusiasm for optimal 
signal processing is somewhat dampened when we come to practical realization. The 
receiver structures for regions 2 and 3 require an operation of the form: sample, 
square-law, and integration. This involved operation cannot be realized easily in an 
optical instrument as, for example, the simple integration of the mutual coherence by 
means of a long exposure photograph can be. Thus it appears that optimal processing, 
at least at present, requires additional off-line computing hardware. If, then, we are 
going to insist on recreating a full picture of the object, the sheer number of spatial 
frequencies to be estimated is very discouraging. On the other hand, high-fidelity pic- 
tures, though aesthetically pleasing, may often be a very expensive method of obtaining 
the actually desired information. For example, if all we really want to know are the 
general dimensions of the object, we can usually obtain these answers simply from an 
estimate of the roll-off characteristics of the magnitude of the spatial Fourier trans- 
form. Thus one of the fundamental engineering tasks is to distinctly define the esti- 
mating assignment. After clearly understanding what is the desired information and 
what in turn are the significant parameters of the spatial Fourier transform that must 
be determined, we can then turn to the hardware aspects of the problem. 

Most likely we shall be unwilling to go to the expense of exploiting all of the available 
-channel diversity, and we shall accept some sacrifice in performance for receiver cost 
reduction. The optical frequency diversity would probably first be discarded to avoid 
the expense of parallel sets of nonabsorption bandpass filters. Likewise we might be 
tempted to discard spatial diversity because each independent spatial sample requires 
a distinct detector element. Time diversity is likely to be the most attractive to exploit. 

The choice of detecting element is also a prime consideration. Photography has 
been an attractive medium for image recording, primarily because of the large spatial 
information capacity. For our applications, however; a direct photoelectric detector 
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such as a photomultiplier is probably more suitable, because of the increased quantum 
efficiency and dynamic range. 

Because of the number of correlation products and averages that must be computed, 
and the need to avoid any systematic biases, it appears that off-line digital processing 
would not only be cheaper but more reliable than analogue schemes. 

5. 3 SUGGESTIONS FOR FURTHER RESEARCH 

The basic aspects of the imaging problem are now sufficiently well developed that 
the emphasis should next be placed on applying the techniques of signal processing to 
problems of actual practical interest. This will require convincing experimental dem- 
onstrations of the merits'of more sophisticated processing over existing techniques. 

Some of the desirable experimental research is possible over the Museum of Science - 
Harvard Observatory path. As we have mentioned, a mercury arc lamp with various 
spatial frequency gratings was used to some extent, but our work with this was primarily 
heuristic. It would be very worth while to continue an experimental verification of our 
results with the use of such a system. One should, however, replace the present lamp 
with a direct-current model in order to eliminate the annoying 120-Hz fluctuations. 

In terms of channel measurements there is, of course, much more that can be done. 
Even the lowly phase -structure measurements should be examined further, probably 
with smaller separations than the minimum of 2 inches that we employed so that the 
saturation characteristics might be better understood. Undoubtedly, the disappearance 
of the 5/3 behavior in our measurements will be disconcerting to some, and it would 
be desirable to obtain confirmation (or perhaps denial) by an independent observer. In 
conducting additional phase structure function measurements, we suggest that a better 
fringe scanner be used so that measurements can also be taken at times of more rapid 
amplitude fading than we could study. 

Other important aspects of the problem which also deserve study, both experimen- 
tally and theoretically, are the angular and optical frequency dependencies of the atmo- 
spheric disturbances. A good physical understanding of just when the channel is 
underspread or overspread would be valuable. Indeed the whole problem of the over- 
spread channel, that is, when the angular correlation scale is too small to be resolved 
by a lens matched to a coherence area, has been left untouched. 

Recently developed lasers operating simultaneously on several wavelengths would 
be very useful tools for studying the optical frequency correlation properties of the dis- 
turbances. 

Our analysis of the problem was primarily classical, and as is the case when 
tackling a problem that is basically quantum in nature, we run into the difficulty that 
strict system optimization is somewhat of a cut-and-try proposition, since, as observ- 
ers, we inevitably affect the experiment* A strict quantum-theory approach might prove 
quite worthwhile, especially if it could indicate directly the best measurement operator 
to employ* 
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Finally, if all of these suggestions are exhausted, the problem of making the pro- 
cessors more practically attractive will probably always be open. It has been twenty 
years since Claude E. Shannon showed the tremendous merits of channel coding, and 
still the problems of efficient decoder design persist. 
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Appendix A 
(Section H) 


A. 1 ERROR BETWEEN SAMPLE AND ENSEMBLE AVERAGES 
OF MUTUAL COHERENCE 


The mean- square error in the sample average mutual coherence is 

M 2 = |r v (x,x’,t)-r v (x ? x*,t)| 2 

= |r v (x,x>,t) | 2 - |r v (x,x\t) | 2 . 


(A. 1) 


The first term is 

I r v <x, x' , t) [ 2 = Ax 4 J dt J dt' E v {x, t) E v (x> , t) E v (x, F) E y (x> , F) . (A. 2) 

AT AT 

The fourth-order expectation can be reduced, since the process E v (x,t) is Gaussian, 
and we are therefore able to rewrite (A. 2) 

| r v <x, x> , t) | : 2 = R l + R 2 + R 3 , (A. 3) 

where 

R 1 = Ax 4 j* dt J dt' E v (x,t) E v (x',t) E v (x,t') E y (x',t') 

AT AT 

R 2 = Ax 4 j* dt df E v (x,t) E v (x',f) E v (x',t) E v (x,f) 

AT AT 

R 3 = Ax j dt j dt' E v (x,t) E v {x,t') E v (x',t) E y <x',F). {A. 4) 

AT AT 

i A . 2 

The first term in expression (A. 4) is simply | F^x.x' ,t)[ .which conveniently cancels 
with the second term of (A. 1). The middle term has zero expectation, because of the 
uniform distribution of the phase angle of each radiation process. Only the final term 
remains,. 

We now revert to the Fourier series 
M(q+1)-1 

E (x,t) = — Y E (x) e - 3 2irmt /AT nAT <t «{n+l)AT (A. 5) 

v 1 / o l—J *** "" 

AT 17 * m=Mq 
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and exploit the following properties: 

1* different optical frequency components are uncorr elated 

2. e J^^t/AT and e +j2inm t/AT ape Qj^hogonal over the AT interval when m =£ 
to find 


m T 


M(q+l)~l 

I l E n,m ®! 2 


(A. 6) 


m=Mq 


Upon setting M = AT Av and rewriting |E (x) | in terms of the object and noise power, 
we obtain 


I e j 2 = AT Av Ax^ 


2a (x,t) 

® v (Q)e v +N v £2 


_ 2a (x\t) 

© (0)e v +N Q 
v' V 


(A. 7) 
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APPENDIX B 


(Section III) 

B. 1 ASYMPTOTIC EXPANSION OF C-R INEQUALITY FOR 
SMALL VALUES OF A 


We expand the right-hand side of Eq. 122 in a power series in A and retain only the 
lowest order terms. The condition of probability p(r/A, cf>) is written 


i /“V [lr| 2 +Ae 2a -2|r| A 1 / 2 cos (b+^-Q)]' 
p(r/A, <j>) = -£7 ( e 


/a, b 


(B. 1) 


where r = |r | e JW and ( ) ^ denotes the averaging over the atmospheric functions a i 

and b^. Setting A 1 / 2 = x, we have 


4A 


Variance A = 


f -i2 9 

NE|^lnp(r/x,<j>) 


where 


0 

eJ In P(r/x, <f>) = 


<^|r j e a cos (b+<j>-0) e x|r| e * cos (b+<f>-0) -x 2 e 2a /2^ 


x |r | e a cos (b+<j>-0) -x 2 e 2a /2 


a. b 


(B. 2) 


(B. 3) 


We now expand expression (B. 3) as a ratio of power series in x as 

3 1n . , .. Y q + V^+ y 2 x 2 ... 

-g^-ln p(r/x, 4) = 5 


1 + 6jX + 6 2 x . . . 


(B.4) 


and find for the numerator coefficients 


Y 0 = < M e a cos(b+<j>-0)) a b 
= j r | D^ Bj^ cos (<{>— 0) 

Vl = ( | r | 2 e 2a cos 2 (b+<j>-0) - e 2a ) &j b 

jr| 2 D 2 jr| 2 D 2 B 2 cos (2<}>-20) - D 2 
= 2 + 2 
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Y 2 = Jr* ) e 3a cos (b+^-Q) + y |r| 2 e 3a cos 3 (b+^-e)^. 

= ^T |r|D 3 B 1 +| |r| 2 D 3 <cos 3 (b+<j>-e)) b 
and for the denominator 

6 X = < |r| e a cos (b+<j>-0) b 

= [r | Dj Bj cos (<t>— 9) 

6 2 = <(xj r | 2e2acos2 ( b+< f >-0 } b 


|r| 2 D 


-D, 


4-^ [1 + b 2 cos ( 2 4>-26)] 


Dividing and squaring the two series gives the representation 

r 0 ~ ^ 

[te p(r/x ' +> J = v 0 + x2 v 0 <v 1 -y 0 y 1 ) 

+ x 2 [2 7o (y 2 +6^v 0 -5 2 Y 0 ~6 1 Y 1 ) ~ (Y^y/] 
+ higher order terms in x. 


<B. 5) 


(B. 6) 


(B. 7) 


This expression must now be averaged over the random variable r which is itself 
dependent upon x* If we express the resultant average as the power series 


Er 


8 lnp(r/x, <j>) 


9x 


= g Q + gjX + g 2 x + . . . , 


2 . 


' (B.8) 


we find that g , the zero-order term in v _ is 
& o T o 


0 °zero-order ^1 ®l* 


(B. 9) 


Since y Q contains only zero- and second-order terms, the coefficient must be 


the zero-order component of 2 Y 0 (Yj“Y 0 6 ^)> where 


2 Y 0 (Y 1 “Y 0 6 1 ) = 2 r Dj B } cos {((>-0) 


r I 2 D_ Irl D„B„ 

— + _£_£ cos (2<j,-28) -D 2 


3 3 3 3 

2 Jr* j D^Bj cos {4>-0). 


(B. 10) 
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2 

Upon Inspection of this expression, we note that g^ must be of order B 1 or less, and 

hence the term g x will only be significant in comparison with g (which also is of order 
2 1 ° 

Bj) when x is of the order of one or more* If we restrict ourselves only to values of x 

much less than one, the first-order term is then negligible. 

The next coefficient in the expansion g^ is composed of 

2 

1. the second-order component of y^ 

2 . the first-order component of 2 Y 0 (Y 1 ”Y 0 5 j) 

3. the zero -order component of 2 Y c (y 2 +5 i V 0 “ S 2^o _S l^l ) + 

2 

All of these components will be of the order of or less, except for part of the zero- 

2 ^ 

order component of (Yj - Y 0 6 i) • Since Bj is very small in cases of practical interest 
(variance of b is very large) , the only significant component of g 2 is g 2 k D 2 . 

An asymptotic approximation to the left-hand side of the C-R inequality can now be 
constructed by using only the first two significant terms g + g 9 x , and yields 

O b 


4A 


Variance A 


n[dX + d^a] 

2 

If we normalize by dividing by A we obtain 

4 


Variance A 


N D 2 B 2 A + D 2 A 2 J 


(B. 11) 


(B. 12) 


B. 2 APPROXIMATION TO C-R INEQUALITY WHEN AMPLITUDE 
IS CONSTANT AND PHASE IS UNIFORM 


Under the conditions of constant amplitude and uniform phase, the conditional proba- 
bility of a particular r. is 


P(r/A, <j>) =■ 


— -|[| r | 2 + AD> /2 ] A l/Z 


Zn 


j_ r * 

2 * J 0 


D^ 2 | r I cos (b+<j>-0) 


dB. 


(B. 13) 


The inner integral can be expressed in terms of a zero-order modified Bessel function 

22 

of the first kind I Q (x) which is independent of <j> and 0 (see Wozencraft and Jacobs ), 
and we obtain 


P( r /A, <j>) 


1 z + ae* 2 ] 

2 6 



(B. 14) 


Differentiating the logarithm of expression (B.14) twice with respect to A gives 
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(B. 15) 


T Ji T 

a 2 o nA 2 A o 

2 " In p(r/A,<j>) = 

BA I 


The function I 0 ( x ) has the power -series expansion 


I o (X) " 1 + 2 2+ 2 2 4 2+ 2 Z 4 Z 6 Z+ "" 


(B. 16) 


From a convexity argument we can show that the first term of the numerator of 
Eq. B. 15 is always positive, and hence 


^-2 lnp(r/A,4,)< f 


Upon letting x = a}' Z 'D^ Z t, Eq. B. 17 becomes 


(B. 17) 


in P<r/A, <j>) 


D 2 M 2 r d i„h 2 


4A dx 


(B. 18) 


Furthermore, it can be shown by using the power-series expansion for I q (x) that 




(B. 19) 


hence. 


2 X) 

P(r/A,<j>) « — Irl 4 . 


(B. 20) 


Averaging over r yields 


a 2 D^A 2 + 8D 2 A + 8D 2 

In p(r/A <j>) 

9A^ ib 


(B. 21) 


and the lower bound to the variance of A then is 


Variance A ~ ^ 

D^A + 8D|A + 8D* 


{B. 22) 


and dividing by A , we finally obtain 
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(B. 23)j. 


Variance A 


N 8 D 2 A 2 + 8 D 3 A 3 + D 4 A 4 * 


B. 3 ASYMPTOTIC EXPANSION FOR C-R BOUND TO PHASE VARIANCE 
IN THE LIMIT OF VERY SMALL SIGNAL LEVEL 

The C-R inequality requires 


In p(r/A, <S>) 

8<j> 


where A 1 / 2 = x, and r = [r[ e JU *. To first order in-x, this behaves simply as 

8 lnp(r/A,«|>) ~<-x|r| e a sin (b+<j>-0) > b 

8cj> * 

~-x|r|D^B^ sin (<fj — 0 ) . (. 

Squaring, averaging over r, and retaining only the lowest order term in x yields 

9 In p(r/A,<f >) 2 ~ x 2 D 2 B^. (: 

d<l> 1 1 

We thus obtain in the limit of very small A 

a 1 

Variance <f> 3 5 — 5 — . (. 

nd^bJa 

B. 4 APPROXIMATION TO C-R INEQUALITY FOR NOISY PHASE 
REFERENCE AND CONSTANT AMPLITUDE 


' 1 r a . , Q . x ! r ! e a cos (b+A-6) -x 2 e 2a /2 

-x r e sm (b+<j>-0) e ' T ' • 


x | r | e a cos (b+«f>— 0 ) -x 2 e 2a /2 


The joint probability density of a pair of samples r. and v. (under the assumption 

ai l /2 ^ ^ 

that the fading coefficient e is replaced by D,/ ) is 


. - 5 |lrl Z +|v| 2 + 2 AD,] / I cos (b+*-i 

(r,v/A, 4 .) = i e Z L 2 J/e 21 


9) + I v | cos (b-y)] 


r = r e J and v- v e J 


The term that must be averaged on b can be rewritten after some trigonometric 
gymnastics as 
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A 1 / 2 D^/ 2 Mr! 2 + ivl 2 + 2|r| |vl cos (<j>+Y - 0)l cos {b+s\ 

Ob = < e /v 


(B. 30) 


The exact functional form of the angle 6 is unimportant, since the average over b is 

independent of 6 anyway. This average is again expressed in terms of the zero-order 

modified Bessel function I (x) as 

o' 


< >b = h i o a1/2d 2 /2 [I 1 *I 2 + M z+2 M M cos(<i>+-e)] / ]. 


It then follows that 


9(|> 


In p(r, v/A,4>) 


72 

d/dy i D (y)" 

2 

V 

2 

3x 


yy) 


ax 

u J 


84> 


where y = x lj ^ 2 , and x » AD.[ | r 1 2 + j v | 3 +2|r | |v | cos (<j>+ -0)]. 


With the aid of inequality (B. 19] we obtain 

12 


— In p(r,v/A,4>) 


= -jA 2 D 2 | r | 2 | v j 2 sin 2 (<j>+-0). 


This is now averaged over r> and v to yield 


— In p(r,v/A,<f>) 
0<|> 


~]2 A 4 D 4 A 3 D 3 A 2 D 2 

+ - + o < 


8 


The final result then is 


(B. 31) 


(B. 32) 


(B. 33) 


(B. 34) 


Variance <j> 


N A 2 +4D^ A 3 +D^ A 4 ^ 


(B. 35) 
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